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EXECUTIVE  SUMMARY 


Water  for  Life:  Alberta’ s Strategy  for  Sustainability  (Government  of  Alberta  2003) 
requires  regular  monitoring,  évaluation,  and  reporting  of  aquatic  ecosystem  health. 
Although  comprehensive  water  quality  information  is  generated  annually  as  part  of 
Alberta  Environment’ s Long-term  River  Network  (LTRN),  comparable,  current 
information  on  sédiment  quality  and  biological  communities  is  not  available. 

In  2006  benthic  invertebrates  and  epilithic  algae  were  collected  from  the  Bow  River  in 
order  to: 

1 . Fill  knowledge  and  information  gaps  identified  by  North  South  Consulting  Inc.  et 
al.  (2007)  with  current  data  for  Bow  River  LTRN  sites. 

2.  Compare  various  benthic  invertebrate  sampling  techniques  for  long-term  river 
monitoring  sites  in  Alberta  and  provide  recommendations  for  future  monitoring  at 
these  sites. 

3.  Evaluate  temporal  changes  and  longitudinal  trends  in  benthic  invertebrate  and 
epilithic  algal  communities,  particularly  as  they  relate  to  nutrient  management  in 
the  Bow  River. 

After  comparing  the  relative  merits  of  data  collected  with  a Neill  cylinder  sampler  and 
kick  nets  of  varying  mesh  sizes,  the  former  was  selected  for  further  monitoring  at  LTRN 
sites.  This  decision  was  based  on  several  considérations.  Maintaining  continuity  with 
Alberta  Environment’ s historical  database  was  viewed  as  important.  Neill  cylinder 
samples  are  quantitative  and  replicated,  and  hence  lend  themselves  easily  to  statistical 
différence  testing  of  population  densities.  Furthermore,  Neill  cylinder  samples  reflect  a 
greater  taxonomie  diversity  than  kick  samples.  Nevertheless,  kick  samples  could  be 
appropriate  in  other  contexts  even  though  they  only  provide  qualitative  information. 

Benthic  invertebrate  data  from  2006  show  evidence  of  perturbation  at  Cochrane,  most 
likely  due  to  sharp  diurnal  fluctuations  in  discharge  (and  stage)  caused  by  the  operation 
of  an  upstream  hydro-electric  dam.  As  well,  Didymosphenia  geminata,  a colonial  diatom 
which  forms  dense  cotton-like  mats  on  the  rocks  contributes  further  to  the  inhospitable 
nature  of  the  physical  habitat  at  this  site.  Zoobenthic  and  epilithic  algal  communities 
provide  clear  evidence  of  nutrient  enrichment  at  Stier’s  Ranch  and  Carseland  downstream 
of  Calgary.  Algal  and  invertebrate  abundances  are  lower  at  Cluny,  but  a diverse 
community  including  sensitive  taxa  is  supported.  The  fauna  at  Ronalane,  near  the  mouth 
of  the  Bow  River,  is  characteristic  of  a depositional  area  with  a sédiment  bottom, 
macrophytes,  and  potentially  lower  dissolved  oxygen  levels.  Longitudinal  changes  in 
benthic  invertebrate  communities  are  similar  in  erosional  and  depositional  areas,  although 
signs  of  stress,  such  as  lower  diversity  indices,  are  apparent  in  depositional  areas  of 
Stier’s  Ranch  and  Carseland. 

There  is  a gap  of  thirteen  years  between  the  current  data  and  the  last  zoobenthic  samples 
collected  from  the  Bow  River.  Based  on  this  limited  record,  it  appears  that  while  the  City 
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of  Calgary’s  nutrient  removal  program  led  to  some  improvements  reflected  in  the  river’ s 
benthic  invertebrate  and  epilithic  algal  communities,  the  degree  of  impairment  may  again 
be  increasing.  These  conclusions  may  be  biased  by  differing  annual  flow  régimes;  in  the 
future,  a more  consistent  sampling  schedule  should  be  adhered  to  so  that  such  variation 
can  be  better  accounted  for. 

Evaluating  these  biomonitoring  data  along  with  sédiment  and  water  quality  chemistry 
will  help  to  further  define  current  conditions  and  past  changes  in  the  Bow  River.  The 
development  of  indicators  and  perhaps  even  a prédictive  model  for  aquatic  ecosystem 
health  should  follow.  Again,  this  work  would  be  of  most  value  if  supported  by  a regular 
and  consistent  sampling  schedule. 
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1.0  INTRODUCTION 


Regular  monitoring,  évaluation,  and  reporting  of  aquatic  ecosystem  health  are 
requirements  highlighted  under  Water  for  Life,  Alberta’ s Strategy  for  Sustainability 
(Government  of  Alberta  2003).  A recent  report  by  North  South  Consulting  Inc.  et  al. 
(2007)  reviewed  available  information  on  aquatic  ecosystem  health  in  Alberta.  Although 
comprehensive  water  quality  information  is  generated  on  a regular  basis  as  part  of 
reporting  on  Alberta  Environment’ s long-term  river  network  (LTRN),  data  for  sédiment 
quality  and  biological  communities  at  the  network  sites  was  not  available. 

The  use  of  benthic  invertebrate  and  epilithic  algal  communities  in  biomonitoring 
programs  is  widespread  and  well  documented  ( e.g .,  Davis  and  Simon  1995,  Barbour  et  al. 
1999).  Although  recent  biomonitoring  data  are  lacking  in  Alberta,  there  is  a substantial 
historical  database  for  benthic  invertebrates  as  well  as  some  epilithic  algae  data  dating 
back  to  the  1980s  (e.g.,  Charlton  1986,  Charlton  et  al.  1986,  Anderson  1991,  Sosiak 
2002). 

Alberta  Environment  established  a long-term  zoobenthos  monitoring  network  in  the 
1970s  and  80s  at  about  20  sites  on  major  provincial  ri  vers.  The  network  was  designed 
with  spring  and  fall  collections  over  5 consecutive  years  alternating  with  5 years  with  no 
collections.  The  expectation  was  that,  over  time,  it  would  be  possible  to  relate  changes  in 
community  composition  to  changes  in  aquatic  ecosystem  quality  (Anderson  1991). 
Sampling  was  not  resumed  as  scheduled  in  1993  because  of  fiscal  constraints  and 
evolving  provincial  monitoring  priorities.  However,  benthic  invertebrate  monitoring 
continued  in  some  rivers  as  part  of  other  projects  (e.g.,  longitudinal  surveys  in  the  Bow 
River  in  1993,  folio wing  some  major  upgrades  to  municipal  wastewater  treatment  in 
Calgary). 

The  Bow  River  is  an  interesting  System  in  which  to  study  responses  of  biological 
communities  to  improvements  in  municipal  wastewater  treatment.  Phosphorus  removal 
technology  was  implemented  at  Calgary’ s two  wastewater  treatment  plants  in  1982-83. 
This  was  followed  by  improvements  to  nitrogen  removal  between  1987  and  1990. 
However,  between  1983  and  2006,  Calgary’ s metropolitan  area  population  grew  from 
620,692  to  991,759  (City  of  Calgary  2007)  and  is  now  over  1,000,000.  Despite 
réductions  in  nutrient  concentrations  due  to  improved  municipal  treatment,  increasing 
population  density  will  lead  to  increased  wastewater  production  which  could  increase 
nutrient  loadings  to  the  Bow  River. 

In  the  Bow  River,  chlorophylle  and  taxonomie  composition  of  epilithic  algal 
communities  as  well  as  biomass  and  macrophytes  species  composition  hâve  been 
monitored  regularly  since  the  early  1980s.  These  data  were  used  to  document  biological 
responses  to  the  management  of  municipal  effluent  loading  from  Calgary  (Charlton  et  al. 
1986,  Sosiak  2002).  According  to  Sosiak  (2002),  changes  to  nutrient  loading  from  the 
City  of  Calgary  resulted  in  a réduction  of  macrophyte  biomass  but  did  not  significantly 
influence  epilithic  algal  growth  in  the  period  studied  (up  to  1996).  Both  the  continued 
phosphorus  removal  (alum)  and  the  scouring  action  of  the  2005  flood  hâve  resulted  in  an 
aquatic  plant  community  that  is  mostly  periphyton  and  which  now  may  be  responding  to 
much  lower  levels  of  phosphorus  (Al  Sosiak,  pers.  communication,  2008). 
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Benthic  invertebrate  monitoring  was  part  of  surveys  conducted  in  1989  and  1992.  Initial 
benthic  invertebrate  responses  to  phosphorus  removal  from  Calgary’s  municipal 
discharges  were  documented  in  Anderson  (1991),  but  benthic  data  hâve  not  been 
evaluated  since. 

In  2006  a pilot  sampling  program  for  both  benthic  invertebrates  and  epilithic  algae  was 
implemented  on  the  Bow  River.  Objectives  of  the  study  were  to: 

1 . Start  filling  knowledge  and  information  gaps  identified  by  North  South 
Consulting  Inc.  et  al.  (2006)  with  current  data  for  Bow  River  LTRN  sites. 

2.  Compare  the  relative  merits  of  various  benthic  invertebrate  sampling  techniques 
for  long-term  river  monitoring  sites  in  Alberta  and  provide  recommendations  for 
future  monitoring  at  these  sites.  Although  Neill  cylinder  samples  hâve  been  the 
standard  sampling  technique  by  AENV,  rapid  assessment  sampling  techniques  are 
gaining  popularity  and  needed  investigation. 

3.  Evaluate  temporal  changes  and  longitudinal  trends  (both  current  and  historical)  in 
benthic  invertebrate  and  epilithic  algal  communities,  particularly  as  they  relate  to 
nutrient  management  in  the  Bow  River. 
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2.0 


METHODS 


2.1  Sampling  Program 

2.1.1  Introduction 

Alberta  Environment  (AENV)  collected  benthic  invertebrate  and  epilithic  algal 
community  samples  from  five  long-term  monitoring  sites  on  the  Bow  River  (Cochrane, 
Stier’s  Ranch,  Carseland,  Cluny,  and  Ronalane  - Figure  1)  in  October  2006.  Method 
comparison  was  one  of  the  objectives  of  this  project;  therefore,  care  was  taken  to 
alternate  the  collection  of  sample  types  to  account  for  habitat  variability  and  to  augment 
the  comparability  of  data  sets.  For  example,  in  erosional  areas,  six  rocks  for  epilithic 
algal  sampling  would  be  collected  first,  followed  by  one  Neill  cylinder  sample,  a ‘90 
second’  kick  sample  with  a 210|im  net,  and  finally  a ’90  second’  kick  sample  with  a 
400pm  mesh.  This  cycle  would  be  repeated  until  ail  samples  were  collected.  Sampling 
proceeded  in  an  upstream  direction,  taking  care  not  to  disturb  areas  yet  to  be  sampled. 

In  addition  to  biological  sampling  described  below,  flow  velocity,  depth,  and  a visual 
assessment  of  substrate  type  were  recorded  from  the  reach  sampled  at  each  site 
(Appendix  1). 

2.1.2  Benthic  Invertebrates 

Alberta  Environment  has  used  a Neill  cylinder  to  collect  lotie  benthic  invertebrates  for 
over  30  years  ( e.g Anderson  1991).  This  method  is  quantitative  as  the  cylinder  covers  a 
known  area  (0.  lm2)  which  allows  estimâtes  of  invertebrate  population  density.  The 
cylinder  is  inserted  securely  in  the  substrate  so  that  a seal  is  established.  Rocks  enclosed 
in  the  cylinder  are  cleaned  by  hand  and  the  substrate  is  then  agitated  with  a pointed 
shovel  until  water  flows  clear  into  the  collecting  net.  Débris  and  invertebrates  are 
washed  into  a net  (mesh  aperture  of  0.210  pm)  attached  to  the  downstream  portion  of  the 
cylinder,  and  from  there  into  a collection  bottle.  In  accordance  with  AENV  protocol 
(AENV  2006),  five  replicate  samples  are  collected  at  each  site. 

Semi-quantitative  kick  net  sampling  methods  are  favoured  by  many  organisations  (United 
States  Geological  Survey,  USGS,  Cuffney  et  al.  1993;  United  States  Environmental 
Protection  Agency,  USEPA,  Barbour  et  al.  1999;  Ontario  Benthic  Biomonitoring 
Network  OBBN,  Jones  et  al  2005).  These  methods  are  primarily  used  in  conjunction 
with  “rapid  assessment”  and  “reference  condition”  approaches  which  usually  involve 
only  one  sample  per  site.  However,  the  reference  condition  approach  requires  individual 
samples  from  numerous  sites  to  provide  sufficient  sample  size  to  define  the 
characteristics  of  reference  sites.  These  samples  may  be  standardized  based  on  sample 
area,  and/or  the  length  of  time  spent  agitating  the  substrate  upstream  of  the  net.  The 
mesh  size  of  the  net  tends  to  vary.  USEPA  (Barbour  et  al.  1999)  and  OBBN  (Jones  et  al. 
2005)  protocols  specify  500pm  openings;  Waterwatch  Australia  (2004)  and  others 
(Taylor  1997)  recommend  a 250  pm  mesh  size. 

Three  methods  were  used  to  collect  benthic  invertebrates  from  erosional  areas  in  the  Bow 
River: 
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1 . Neill  cylinders  with  a 210|im  mesh  size  and  a 0. 1 m2  sampling  area  Five  (5) 
replicate  samples  per  site  as  per  Alberta  Environment  (2006). 

2.  Kick  samples  taken  with  a 400pm  mesh  size,  total  time  of  3 minutes.  A single 
sample  per  site  with  two  extra  samples  at  2 sites  for  QA/QC  purposes. 

3.  Kick  samples  taken  with  a 210pm  mesh  size,  total  time  of  3 minutes.  A single 
sample  per  site  with  two  extra  samples  at  2 sites  for  QA/QC  purposes. 

Each  kick  sample  was  usually  made  up  of  two  subsamples  timed  at  90  seconds  and 
pooled  together  to  form  a composite.  The  net  was  emptied  more  often  if  fines  were 
clogging  it,  but  the  total  sampling  time  was  always  3 minutes.  Ail  invertebrate  samples 
were  preserved  shortly  after  collection  with  4%  formaldéhyde. 

In  depositional  areas  a kick  sample  taken  with  a 210pm  mesh  size  was  collected  over  90 
sec. 

2. 1.3  Epilithic  Algae 

Epilithic  algae  samples  for  chlorophyll-a  analysis,  as  well  as  samples  for  taxonomie 
analysis  were  collected  from  erosional  areas.  Epilithic  algae  for  chlorophyll-a 
détermination  were  scraped  from  rocks  using  the  template  method  (Alberta  Environment 
2006).  Scrapings  from  a 4 cm2  template  were  taken  from  each  of  3 rocks  to  form  a 
composite  sample.  Algal  material  was  placed  on  a GF/C  filter,  sprinkled  with  MgCC>3, 
and  then  wrapped  in  aluminum  foil  and  frozen.  Three  replicate  composite  samples  were 
collected  at  each  site,  and  up  to  6 additional  samples  were  taken  at  some  sites  for  QA/QC. 

Epilithic  algae  for  taxonomie  analysis  were  also  obtained  using  the  template  method,  but 
in  this  case  scrapings  from  9 rocks  (4  cm"  scraping  each)  were  combined  to  form  one 
composite  sample.  The  sample  was  preserved  shortly  after  collection  with  Lugol’s 
solution  and  5 drops  of  formaldéhyde. 


2.2  Laboratory  Analyses 

2.2.1  Benthic  Invertebrates 


The  zoobenthic  samples  were  stained  with  Rose  Bengal  (Mason  and  Yevich  1967)  then 
sorted  through  three  sieves  with  mesh  apertures  of  2,  1,  and  0.2  mm.  The  coarse  fraction 
and  fine  fractions  were  sorted  under  a dissecting  microscope  (magnification  range  6 to 
50X).  It  was  necessary  to  sub-sample  the  fine  fractions  (residue  on  the  1 and  0.2  mm 
screens)  of  some  samples  that  contained  large  numbers  of  organisms.  Sub-sampling  was 
performed  using  the  Imhoff  cône  method  described  by  Wrona  et  al.  (1982). 

Specimens  were  identified  to  genus  or  species  where  possible,  according  to  Edmunds  et 
al.  (1976),  Baumann  et  al.  (1977),  Wiggins  (1977),  Pennak  (1978),  Merritt  and  Cummins 
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(1984,  1996),  Clifford  (1991),  Thorp  and  Covich  (2001),  and  others  using  the  most 
current  taxonomie  désignations  available. 

Raw  benthic  invertebrate  data  are  presented  in  Appendix  2 


2.2.2  Epilithic  Algae 

Chlorophyll-tf  and  phaeophytin-^,  a dégradation  product  of  chlorophyll,  were  determined 
fluorometrically  after  acetone  extraction  at  Alberta  Research  Council,  Vegreville. 

Results  are  reported  as  mg/m2. 

Non-diatoms  (soft  algae)  and  diatoms  were  analyzed  separately.  Depending  on  their 
concentration,  non-diatom  samples  were  diluted  first.  To  détermine  the  appropriate 
dilution,  the  original  samples  were  screened  to  assess  the  densities  of  algae  and  non-algal 
matter  (débris  and  particulate  matter).  Aliquots  of  the  appropriately  diluted  samples  were 
allowed  to  settle  overnight  in  sédimentation  chambers  following  Utermôhrs  procedure 
described  in  Lund  et  al.  (1958).  Algal  units  were  counted  from  a minimum  of  four 
transects  on  a Zeiss  Axiovert  40  CFL  inverted  microscope.  Counting  units  were 
individual  cells,  filaments,  or  colonies  depending  on  the  organization  of  the  algae.  Both 
diatoms  and  non-diatoms  were  counted.  For  soft  algae,  between  250  and  300  units  were 
counted  at  500X  magnification;  a number  of  transects  were  scanned  at  250X  for  larger 
algae.  For  diatoms,  a minimum  of  250  was  set  as  the  target.  At  this  stage,  diatoms  were 
not  identified  to  species  or  genus,  but  recorded  as  “diatoms”,  and  were  later  identified  to 
species  from  prepared  slides. 

Préparation  of  diatom  slides  consisted  of  digesting  sub-samples  using  concentrated  nitric 
acid  and  hydrogen  peroxide  and  washing  several  times  (by  centrifuging)  with  distilled 
water.  A few  drops  of  the  diatom  slurry  were  placed  on  a cover  slip  and  allowed  to 
evaporate  overnight.  Once  dry,  the  diatoms  were  mounted  in  Naphrax  and  identified 
using  1000  to  1500  X magnifie ations  (under  oil  immersion)  on  a Zeiss  Axioskop  40 
compound  microscope.  A minimum  of  500  diatom  frustules  were  counted  on  each  slide. 
The  diatom  counts  on  the  slides  were  converted  to  density  based  on  the  number  of 
transects  covered  during  the  fresh  (Utermôhl)  counts. 

Biomass  was  calculated  from  recorded  abundance  and  spécifie  biovolume  estimâtes, 
based  on  géométrie  solids  (Rott  1981),  assuming  unit  spécifie  gravity.  The  bio volume  of 
each  species  was  estimated  from  the  average  dimensions  of  10  to  15  individuals.  The 
biovolume  of  colonial  taxa  was  based  on  the  number  of  individuals  in  a colony.  Ail 
calculations  for  cell  concentration  (units/cm2)  and  biomass  (pg/cm2)  were  performed  with 
Hamilton’s  (1990)  computer  program. 

Taxonomie  identifications  of  soft  algae  were  based  primarily  on  Anton  and  Duthie 
(1981),  Entwisle  et  al,  (2007),  Findlay  and  Kling  (1976),  Huber-Pestalozzi  (1961,  1972, 
1982,  1983),  Tikkanen  (1986),  Prescott  (1982),  Whitford  and  Schumacher  (1984), 
Starmach  (1985),  Komarek  & Anagnostidis  (1998,  2005),  and  Wehr  and  Sheath  (2003). 
Diatom  identifications  were  based  primarily  on  the  following  texts  and  supplemented 
with  other  publications:  Krammer  and  Lange-Bertalot  (1986,  1988,  199 la, b),  Reavie  and 
Smol  (1998),  Cumming  et  al.  (1995),  Bahls  (2004),  Camburn  and  Charles  (2000),  Fallu 
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et  al.  (2000),  Patrick  and  Reimer  (1966,  1975),  Siver  and  Kling  (1997),  and  Siver  et  al 
(2005). 

Raw  epilithic  algal  data  are  presented  in  Appendix  3. 


2.3  Data  Analyses 

2.3. 1 Benthic  Invertebrate  Data 

2.3. 1.1  Metrics 

The  use  of  benthic  invertebrates  as  indicators  of  aquatic  ecosystem  health  is  well 
established.  Characteristics  such  as  feeding  method,  habit,  habitat  preference,  pollution 
tolérance  (primarily  tolérance  to  organic  pollution)  and  environmental  sensitivity  can 
account  for  the  distribution  of  benthic  invertebrates  and  can  be  used  to  illustrate  changes 
due  to  both  natural  and  man-made  factors.  Many  metrics,  based  on  invertebrate  counts 
(individuals  and  taxa),  proportions,  or  more  complex  indices  hâve  been  derived  to 
summarize  invertebrate  data  in  order  assess  river  conditions  (Merritt  and  Cummins  1996, 
Barbour  et  al.  1999,  Mandaville  2002).  Several  such  metrics  were  computed  from  raw 
data  for  this  report.  These  metrics,  along  with  their  rational  and  formulae,  if  appropriate, 
are  listed  here. 

Abundance  - The  total  count  of  individuals  in  a sample. 

Richness  - The  number  of  unique  taxa  in  a sample.  Richness  was  calculated  at  the 
species,  genus,  and  family  level  for  this  report.  A greater  the  number  of  taxa  usually 
suggests  a healthier  aquatic  community. 

Diversity  - Diversity  combines  the  concepts  of  abundance  and  richness.  Two  common 
diversity  indices  were  calculated  for  this  report  (Béais,  Gross,  and  Harrell  1999,  2000): 

Simpson  Diversity  Index  (D)  - 

D Ef.  i »,(»,  - 1 ) 

N (N  - 1) 

. . . where  S is  the  number  of  taxa  (species),  N is  the  number  of  organisms,  and  n is  the 
number  of  organisms  in  the  “ith”  taxon.  This  index  is  usually  presented  as  1-D  so  that  the 
results  are  more  intuitive  (i.e.  the  lower  the  number,  the  lower  the  diversity). 

Shannon  Diversity  Index  (H)- 

h=-2pAüps. 

j= î 

...where  S is  the  number  of  taxa  in  the  community,  and /?,  is  the  proportion  of  individuals 
in  the  “ith”  taxon.  H increases  with  increasing  diversity. 


An  Assessment  of  Benthic  Invertebrates  and  Epilithic  Algae  at  Long-term  Monitoring 
Sites  in  the  Bow  River  (Fall  2006) 


6 


“Number  of  unique  zoobenthos  group s ” - The  number  of  species,  généra  or  families 
within  a particular  group  of  invertebrates,  whether  defined  taxonomically,  by  functional 
feeding  group,  habit/behaviour,  habitat  preference,  or  pollution  tolerance/sensitivity. 
Useful  measures  within  this  category  include: 

Number  ofEPT  taxa  - Members  of  the  orders  Ephemeroptera,  Plecoptera,  and 
Trichoptera  are  generally  considered  to  be  sensitive  to  pollution,  preferring  oxygen- 
rich  habitats.  Richness  in  these  groups  should  be  higher  when  water  quality 
conditions  are  better. 

Number  of  ET O taxa  - This  is  similar  to  the  previous  measure,  but  includes  Odonata, 
instead  of  Plecoptera,  which  are  also  considered  to  be  sensitive. 

Number  ofClinger  taxa  - These  organisms  hold  their  position  on  the  bottom  substrate 
in  flowing  water.  Their  numbers  are  lower  in  the  presence  of  excessive  algal  growth 
or  sédimentation. 

Number  of  Pollution-Intolérant  Taxa  - ‘Pollution’  refers  primarily  to  nutrient  or 
organic  enrichment.  Organisms  were  assigned  a value  between  0 (very  intolérant)  to 
10  (very  tolérant)  based  on  the  concept  of  Hilsenhoff  (1987,  1988)  as  compiled  by 
Mandaville  (2002).  For  this  report,  a conservative  approach  was  used  where  species 
with  rankings  from  0 through  2 were  counted  as  intolérant. 

Proportional  Abundance  - The  number  of  a particular  group  of  invertebrates  in 
comparison  to  the  total  number  of  individuals  in  a sample.  Examples  of  metrics  in  this 
category  include: 

Proportion  of  Chironomidae  - Higher  proportions  of  midge  larvae  are  expected  at 
more  impacted  sites. 

Proportion  ofEPT  - Higher  proportions  of  these  organisms  are  expected  at  less 
impacted  sites. 

Proportion  ofFiltering  Collectors  - These  organisms  will  be  found  in  greater 
proportions  at  sites  with  more  fine  particulate  organic  matter  in  the  water  column. 

Proportion  of  Burrowers  - These  organisms  will  be  found  in  greater  numbers  where 
the  riverbed  is  covered  with  sédiment. 

Proportion  with  Erosional  Preference  - These  invertebrates  should  occur  in  greater 
proportions  in  erosional  (riffle)  areas  (rather  than  in  depositional  or  sédiment 
habitats). 

Percent  Contribution  by  Dominant  Taxon  - A large  prevalence  of  a single  type  of 
organism  in  a sample  suggests  an  impacted  site.  In  this  report,  the  dominant  order  of 
insects  was  compared  to  total  insects. 

Ratio  ofEPT  to  Chironomidae  - An  indication  of  community  balance,  this  ratio  is  lower 
under  conditions  of  environmental  stress. 

Ratio  of  H\drops\chidae  to  Trichoptera  - Hydropsychids  are  considered  to  be  less 
sensitive  than  most  caddisflies. 
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Ratio  of  Scravers  to  Filterers  - This  ratio  compares  two  different  feeding  strategies  to 
reflect  different  environmental  conditions.  Scrapers  prefer  areas  with  lots  of  algae  and 
plant  material,  while  filters  require  fine  particulate  material  in  the  water  column. 

Filterers  may  also  be  used  to  assess  pollutants  adsorbed  to  particulates,  although  they  are 
not  necessarily  sensitive  to  such  pollutants.  This  ratio  did  not  provide  much  information 
for  the  Bow  River  sites. 

Biotic  Index  (BI%  Svecies ) and  Family  Biotic  Index  (FBI)  - (Modified  from  Hilsenhoff 
1987,  1988).  Pollution  tolérance  values  from  0 to  10  were  assigned  at  the  species  or 
family  level  according  to  tables  compiled  by  Mandaville  (2002).  The  indices  were  then 
calculated  according  to  the  following  formula: 

BI  or  FBI  = Z(ni*ti)/(N) 

. . . where  the  number  of  individuals  in  the  “zth”  taxon,  U is  the  tolérance  value  of  the  zth 
taxon,  and  N is  the  total  number  of  organisms  in  the  sample.  A lower  index  value  implies 
better  water  quality. 

Biological  Monitoring  Working  Partx  (BMWP)  - (Modified  for  North  America  by 
Mackie  2004).  This  procedure  assigns  scores,  primarily  at  the  family  level,  between  1 
and  10,  where  1 is  least  sensitive  to  pollution  and  10  indicates  a high  sensitivity.  The 
scores  for  ail  families  are  then  added  together.  A higher  number  indicates  better  aquatic 
conditions.  Scores  for  this  index  were  obtained  from  Mandaville  (2002). 

Average  Score  Per  Taxon  (ASPT)  - The  BMWP  (sensitivity)  for  a sample  is  divided  by 
the  number  of  taxa  used  to  calculate  the  BMWP.  This  éliminâtes  the  effect  of  family 
richness,  yielding  a number  between  1 and  10;  the  higher  the  ASPT,  the  better  the 
conditions. 

2.23.2  Sampling  Method  Comparison 

The  raw  numbers  and  20  metrics  selected  from  Section  2.2.3. 1 were  used  to  compare  the 
three  sampling  methods  based  on  two  broad  criteria:  whether  similar  patterns  were 
exhibited  among  sites  and  whether  the  kick  samples  fell  within  the  variability  exhibited 
by  the  replicate  Neill  cylinder  samples.  The  latter  is  particularly  important  if  valid 
historical  comparisons  are  to  be  made  should  the  Neill  cylinder  method  eventually  be 
replaced  by  kick  samples. 

Method  comparisons  were  accomplished  both  statistically  (11  metrics)  and/or  graphically 
(20  metrics).  Because  five  replicates  are  available  for  the  Neill  cylinder,  the  data  can  be 
treated  in  two  ways,  either  using  the  mean  (or  other  measure  of  central  tendency)  of  a 
metric  calculated  for  each  of  the  five  samples,  or  by  calculating  the  metric  on  the 
averaged  (or  summed)  data.  The  results  can  be  very  different  for  metrics  dealing  with 
counts  (of  species  or  other  taxonomie  levels)  but  are  more  similar  for  proportional  data 
(ratios,  percentages).  The  two  ways  of  dealing  with  the  data,  as  described  above,  are 
designated  as  “Average  Neill”  and  “Total  Neill”  in  the  Figure  2 and  Appendix  4. 
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The  Neill  replicates  were  compared  statistically  to  the  individual  kick  and  total  Neill 
(a.k.a.  Neill  mean  metric)  by  using  the  “bootstrap”  method,  which  estimated  the  mean 
différence  and  95%  confidence  bounds  over  a course  of  1000  itérations  for  each  metric. 
The  bootstrap  offers  a well-established  method  for  estimating  parameters  and  their 
standard  errors  when  larger  samples  are  not  available  (Manly  1991,  2001;  Efron  and 
Tibshirani  1986,  1993). 

The  protocol  was  as  follows: 

1.  Sample  5 values  from  field  data  (replicates),  with  replacement. 

2.  Calculate  the  mean. 

3.  Subtract  the  reference  value  (210  kick,  400  kick,  Neill)  from  the  mean  (#2). 

4.  Store  the  resulting  différence. 

5.  Repeat  the  above  1000  times. 

6.  Calculate  the  upper  and  lower  95%  percentiles  around  the  1000  différences. 

7.  Decision:  do  the  95%  bounds  contain  “0”  (zéro)?  If  so,  this  is  good  evidence 
that  the  respective  sample  method  (210  kick,  400  kick,  or  Neill  Mean)  is  a 
good  fit,  with  95%  confidence,  to  the  Neill  replicates. 

The  total  number  of  factorial  combinations  in  this  exercise  was:  3 methods  x 5 locations 
x 21  metric  = 315. 

Ail  bootstrapping  was  done  using:  Resampling  Stats,  version  5.0.2,  Resampling  Stats, 
Inc.,  Arlington,  VA,  USA,  22201.  The  reasoning  for  performing  1000  itérations  is 
contained  in  Manly  (1991). 


2.2.33  Assessment  of  Longitudinal  and  Temporal  Trends 

In  order  to  compare  the  Bow  River  zoobenthos  data  collected  in  autumn  2006  to  other 
data  also  collected  in  autumn  at  the  same  sites  (i.e.,  1983,  1984,  1985,  1986,  1987,  1989, 
and  1993),  some  changes  to  taxonomie  groupings  had  to  be  made.  Identifications  in  the 
two  previous  décades  were  often  not  as  detailed  as  those  done  for  the  2006  samples; 
occasionally,  identifications  were  at  relatively  high  levels,  family  or  order.  This 
différence  in  level  of  detail  had  the  greatest  impact  on  results  for  the  family 
Chironomidae.  These  organisms  were  almost  always  identified  to  genus  for  the  2006 
samples  but  were  left  at  tribe  or  sub-family  in  the  earlier  samples.  Changes  to  taxonomie 
conventions  in  the  13  years  between  1993  and  2006  also  led  to  différences  between  the 
recent  and  earlier  data  sets.  In  particular,  mayflies  included  historically  in  the  genus 
Baetis  hâve  been  split  into  the  following:  Acerpenna,  Diphentor,  Fallceon,  Procleon, 
Plauditus,  Acentrella,  and  Baetis,  while  the  historical  genus  Psuedocoleon  was  included 
with  Baetis  (W.J.  Anderson,  pers.  comm.  2007,  Merritt  and  Cummins,  2nd  and  3rd 
Editions).  Therefore,  the  data  had  to  be  altered  to  achieve  taxonomie  consistency  among 
years. 

The  raw  data  and  the  derived  metrics  can  be  found  in  Appendix  5. 

Statistical  analyses,  usually  performed  on  transformed  data  ( e.g .,  log(x+l),  arcsin(x)), 
consisted  of  principal  components  analysis  (PCA)  and  step-wise  discriminant  analysis  for 
ail  sites  and  years;  analysis  of  variance  (ANOVA)  for  changes  in  mean  values  over  sites 
for  each  year  ( longitudinal  différences ),  changes  in  mean  values  over  years  for  each  site 
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(, temporal  différences),  and  changes  in  PC  scores  and  discriminant  function  (DF)  scores 
for  the  five  sites  over  time;  and  Bonferroni  pair-wise  comparisons  among  years  and  sites 
for  ail  scores  and  means. 


2.3.2  Epilithic  Al  gai  Data 

Longitudinal  trends  for  chlorophylle  data  and  various  metrics  based  on  the  taxonomie 
data  were  examined  graphically.  Metrics  included  biomass,  density,  number  of  species 
and  généra,  and  % biomass  distributed  over  the  various  algal  groups  for  each  site,  as  well 
as  the  Shannon  diversity  index.  In  addition,  diatom  species  were  classified  as  preferring 
low  or  high  total  phosphorus  concentrations,  and  low  or  high  total  nitrogen 
concentrations,  according  to  the  work  of  Potapova  and  Charles  (2007). 

Pearson  corrélations  were  performed  to  examine  the  relationship  between  various 
invertebrate  and  epilithic  algal  metrics.  A limited  amount  of  historical  epilithic 
community  and  chlorophylle  data  were  also  available  for  comparison. 
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3.0 


RESULTS  AND  DISCUSSION 


3.1  Comparison  of  Benthic  Invertebrate  Sampling  Methods 

3.1.1  Numbers  of  Organisais  and  Taxa 

Because  the  Neill  cylinder  samples  are  quantitative  and  the  kick  samples  are  not,  the 
actual  number  of  organisms  counted  cannot  be  directly  compared.  However,  the  number 
of  organisms  collected  can  be  important  in  that  too  small  a sample  may  not  adequately 
show  the  diversity  of  the  benthic  community,  and  too  large  a sample  présents 
énumération  difficultés.  If  the  collection  methods  are  équivalent,  however,  the 
proportions  of  the  various  groups  of  organisms  collected  should  be  similar,  even  if  the 
numbers  are  not  comparable.  The  Neill  average  had  the  lowest  abundances,  while  the 
total  Neill  and  210  kicks  had  the  highest  abundance.  The  fewest  organisms  were 
collected  at  Cochrane  in  ail  cases.  The  site  with  the  most  organisms  varied  with  method; 
the  largest  number  of  organisms,  73675,  was  collected  in  the  210  kick  at  Stier’s  Ranch 
(Figure  2.1).  Counts  in  the  Neills  were  statistically  different  from  those  in  the 
corresponding  kick  samples  at  ail  sites  (Appendix  6). 

An  important  indicator  of  aquatic  ecosystem  health  is  the  number  of  unique  taxa 
collected.  Samples,  therefore,  must  accurately  reflect  the  variety  of  organisms  présent. 
For  ail  methods,  the  lowest  number  of  taxa  was  collected  at  Cochrane.  The  highest 
number  of  taxa  was  collected  at  Carseland  for  ail  methods  except  for  the  400  kick.  The 
individual  Neill  samples,  which  tend  to  collect  fewer  organisms  than  the  kicks,  had  the 
least  number  of  taxa.  When  the  results  of  the  replicates  were  summed,  however,  the 
Neills  had  the  highest  number  of  taxa  - between  22  and  30  more  per  site  than  the 
individual  samples.  The  second  highest  number  of  taxa  was  usually  collected  by  the  210 
kicks  (Figure  2.2).  The  richness  of  the  kicks  and  of  the  summed  Neills  could  almost 
always  be  considered  to  be  different  than  that  of  the  individual  Neills,  with  a few 
exceptions  at  Stier’s  Ranch  and  Cluny  for  the  210  kicks. 

The  different  collection  methods  produced  similar  longitudinal  patterns  for  many  sub-sets 
of  organisms,  such  as  number  of  chironomid  généra.  The  greatest  numbers  of 
chironomid  généra  were  found  at  Carseland  by  ail  methods.  As  for  overall  taxon 
richness,  the  lowest  numbers  were  from  the  average  Neills  and  highest  were  from  the 
total  Neills  (Figure  2.3).  For  Ephemeroptera,  Plecoptera,  and  Trichoptera  (EPT)  généra, 
ail  methods  found  the  number  of  unique  généra  to  be  highest  at  Cluny.  While  the  total 
Neills  tended  to  hâve  the  highest  counts  overall,  the  400  kick  samples  had  higher  counts 
than  the  210  samples  and  the  average  Neills  had  the  lowest  (Figure  2.4). 

It  is  possible  that  the  behaviour  of  the  benthic  invertebrates  {Le.,  whether  they  burrow, 
climb,  cling,  swim,  etc.)  might  affect  their  capture  by  the  Neill  cylinder  compared  to  kick 
nets.  For  clinger  taxa,  a similar  longitudinal  pattern  was  seen  for  ail  methods,  with  more 
clinger  taxa  being  found  at  Carseland  than  anywhere  else  (Figure  2.5).  For  taxa  classified 
primarily  as  swimmers,  however,  similarities  among  the  methods  were  not  apparent 
(Figure  2.6).  These  results  may  reflect  a sampler-type  bias,  an  overall  challenge 
associated  with  collecting  swimmers,  or  the  fact  that  the  number  of  taxa  in  this  habit 
category  was  low. 
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The  benthic  invertebrates  were  also  categorised  by  functional  feeding  group.  Unlike 
some  of  the  previous  examples,  the  numbers  of  collector-gatherer  taxa  for  both  averaged 
and  total  Neill  samples  were  higher  than  for  the  kick  net  samples.  Lowest  counts  overall 
were  found  at  Cochrane;  for  the  Neills,  the  highest  number  of  gatherers  were  found  at 
Carseland,  while  for  the  kicks,  numbers  were  highest  at  Ronalane  (Figure  2.7). 

Relatively  low  counts  in  the  individual  Neill  samples  occurred  for  the  predator  taxa;  as  in 
the  previous  case,  however,  the  kick  samples  were  much  more  similar  to  each  other  than 
they  were  to  either  the  averaged  or  total  Neills  (Figure  2.8). 

Both  kicks  and  Neills  showed  similar  site-to-site  patterns  for  taxa  classified  by  habitat 
preference  (erosional,  depositional,  etc.).  The  most  taxa  with  an  erosional  preference 
were  found  at  Carseland  and  the  least  were  found  at  Cochrane,  independent  of  the 
method  used.  Overall,  the  greatest  variety  of  erosional  taxa  was  présent  in  the  total  Neill 
samples  (Figure  2.9). 

Ail  sampling  methods  showed  similar  patterns  among  sites  for  taxa  intolérant  to  pollution 
(i.e.,  those  rated  as  0,  1,  or  2 on  a scale  from  0 to  10).  This  occurred  even  though 
relatively  small  numbers  of  these  taxa  were  found  (Figure  2.10).  The  number  of 
intolérant  taxa  in  the  400  kick  samples  was  consistently  higher  than  that  in  the  210  kick 
samples,  despite  the  fact  that  the  210  kicks  collected  more  invertebrates  in  total. 

3. 1.2  Relative  Proportions  of  Organisais 

Calculating  the  proportion  of  various  organisms  to  the  total  collected  or  to  a sub-set  of 
the  collection  should  reduce  some  of  the  différence  among  methods  that  is  due  to  sample 
size.  For  the  Oligochaeta,  longitudinal  patterns  and  percent  contributions  were  similar 
for  the  average  and  total  Neill  samples  and  the  400  kick  samples  at  ail  sites.  At 
Carseland  and  Cluny,  the  percentages  of  the  210  kick  samples  made  up  of  aquatic  worms 
were  double  the  percentages  collected  by  the  other  methods  (Figure  2. 1 1).  This  finding 
may  reflect  the  fact  that  a smaller  mesh  size  is  more  effective  at  collecting  smaller 
organisms;  at  these  sites,  the  oligochaetes  were  almost  exclusively  from  the  family 
Naididae,  which  tend  to  be  quite  small.  The  Neill  and  400  kick  also  tracked  each  other 
more  closely  than  the  210  kick  samples  for  the  percentage  of  Diptera  (including  the 
Chironomidae),  although  the  overall  trend  was  similar  for  ail  samples  (Figure  2.12). 

The  dominance  of  one  type  of  organisms  in  a sample  can  be  an  indication  of  an  impacted, 
or  a less  complex,  benthic  environment.  The  contribution  by  the  dominant  insect  order  at 
each  site  (as  a percentage  of  total  insects)  showed  doser  correspondence  between  the 
Neill  samples  and  the  210  kick  samples  than  the  400  kick  samples  (Figure  2. 13).  The 
percentage  contribution  was  always  lower  in  the  400  kicks  than  the  individual  Neill 
samples,  and  significantly  so  at  four  of  the  five  sites  (Appendix  6).  In  most  samples, 
dipterans  dominated,  primarily  as  a function  of  number  of  Chironomidae.  The  only 
exceptions  were  for  both  Cochrane  kicks,  in  which  Hemiptera  (Corixidae)  were  most 
plentiful  and  for  the  400  kick  at  Cochrane  for  which  the  high  Diptera  counts  were  due  to 
Simuliidae. 

The  ratio  of  EPT  to  Chironomidae  can  be  used  to  indicate  environmental  stress,  with 
higher  ratios  (relatively  more  EPT)  indicating  less  stress.  In  this  case,  the  400  kick 
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samples  gave  much  higher  EPT/Chironomidae  ratios,  particularly  at  Stier’s  Ranch  and 
Cluny  (Figure  2.14).  Not  only  were  fewer  small  and/or  filiform  organisms,  such  as 
chironomids,  collected  by  the  400  mesh  nets,  but  it  appears  that  more  EPT  organisms 
were  collected  as  well  (see  Figure  2.4).  This  finding  was  confirmed  by  the  statistical 
exercise.  With  the  exception  of  Ronalane,  either  one  or  both  of  the  proportion  of  EPT 
and  the  EPT/Chironomidae  ratio  were  significantly  higher  for  the  400  mesh  kick  samples 
than  for  the  individual  Neills  (Appendix  6).  Perhaps  greater  flow  through  the  wider  mesh 
offers  less  résistance  to  larger  and  mobile  insects.  Slower  flow  through  the  finer  mesh, 
slowed  further  by  a greater  chance  of  the  mesh  being  obstructed  by  débris,  may  allow 
more  active  animais  a chance  to  escape  capture. 

The  longitudinal  patterns  for  the  percentage  of  organisms  that  cling  to  the  substrate  were 
similar  for  ail  the  collection  methods,  with  the  percentage  of  clinger  taxa  highest  for  ail 
samples  at  Stier’s  Ranch  and  lowest  at  Ronalane.  At  each  site,  the  percentage  of  clingers 
tended  to  be  lower  for  the  210  kicks  than  the  other  methods  (Figure  2. 15).  As  suggested 
previously,  increased  résistance  to  flow  through  the  210  mesh  net  may  hâve  made  capture 
of  clinger  taxa  slightly  more  difficult.  Swimming  taxa,  the  actual  numbers  of  which  were 
difficult  to  interpret,  clearly  made  up  a larger  percentage  of  ail  samples  taken  at  Cochrane 
than  elsewhere.  The  percentage  of  swimmers  was  much  higher  in  the  210  kick  than  in 
other  samples  at  this  site.  Over  ail  sites,  however,  the  proportion  of  swimmers  was 
slightly  higher  in  the  400  kicks  (Figure  2.16). 

The  percentage  of  predator  taxa  was  higher  for  both  types  of  kick  samples  than  for  the 
Neill  samples  at  Cochrane  and  Stier’s  Ranch.  Sample  results  were  not  as  different  at  the 
three  downstream  sites  (Figure  2.17).  The  ratio  of  scraper  taxa  to  filtering  collectors, 
however,  was  similar  for  ail  sample  types  at  the  first  three  sites.  At  Cluny,  both  kicks 
differed  from  the  Neills  (210  lower,  400  higher)  and  at  Ronalane  the  400  kick  had  a 
higher  scraper  to  filterer  ratio  than  the  Neills  (Appendix  6). 

The  210  kick  samples  at  Cochrane  and  Stier’s  Ranch  appeared  to  be  different  from  both 
the  Neills  and  the  400  kicks  for  taxa  preferring  depositional  habitats.  At  Cochrane,  the 
210  sample  had  the  highest  percentage  of  depositional  taxa  compared  with  the  other 
sampling  methods,  while  at  Stier’s  Ranch  the  210  sample  had  the  lowest  percentage 
(Figure  2.18).  Statistically,  the  210  kick  samples  at  Cluny  and  Ronalane  were  also 
different  than  the  Neills  (Appendix  6).  For  taxa  with  a preference  for  an  erosional 
habitat,  the  proportion  found  in  the  400  kicks  was  different  than  that  of  the  individual 
Neills  at  four  of  the  five  sites  (lower  on  3 of  these  occasions),  while  the  proportion  found 
in  the  210  kicks  was  higher  than  the  Neills  at  Cluny  and  Ronalane,  and  lower  at 
Cochrane. 

3.1.3  Indices 

The  Biotic  Index  is  calculated  by  assigning  pollution  tolérance  values  from  0 (intolérant) 
to  10  (tolérant)  to  ail  taxa  (species  in  this  case)  to  détermine  an  average  tolérance  for  the 
community.  Longitudinal  patterns  were  generally  similar,  with  the  highest  or  second 
highest  scores  for  each  method  occurring  at  Cochrane  and  Ronalane,  and  the  lowest 
scores  (more  favourable  conditions  for  invertebrates)  for  each  method  occurring  at  Stier’s 
Ranch.  However,  the  kick  samples  gave  index  values  that  were  quite  different  from  the 
Neill  samples  and  often  very  different  from  each  other  (Figure  2.19  and  Appendix  6).  It 
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is  apparent  that  the  meaning  of  the  index  values  must  be  calibrated  for  each  method,  and 
need  to  be  reviewed  in  the  context  of  actual  rather  than  assigned  taxa  tolérances,  and 
actual  river  conditions.  For  example,  despite  pronounced  diurnal  fluctuations  and 
critically  low  DO  in  summer  at  Stier’s  Ranch,  the  abundant  food  supply  resulting  from 
enrichment  at  this  site  still  allows  the  establishment  of  large  invertebrate  populations. 

The  Shannon  Index  uses  the  number  of  individuals  in  each  taxon  présent  in  a community 
to  give  an  indication  of  diversity.  A higher  value  represents  higher  biotic  diversity.  The 
index  values  differed  for  each  collection  method,  giving  different  impressions  of  the 
pattern  of  diversity  at  each  site.  The  highest  index  scores  for  the  kick  methods  were 
achieved  at  Cochrane,  while  the  Neill  methods  produced  the  lowest  scores  at  this  site 
(Figure  2.20).  As  for  the  biotic  index,  the  diversity  index  values  should  not  be  compared 
among  different  collection  methods  and  need  to  be  evaluated  in  the  context  of  other  data 
(Appendix  6). 

3.1.4  Other  Considérations 

As  well  as  being  able  to  capture  a représentative  sample  of  benthic  invertebrates,  an  idéal 
collection  method  should  be  easy  to  use  and  be  both  time-  and  cost-effective.  Kick  nets 
and  Neill  cylinders  are  both  fairly  easy  to  use.  The  nets  tend  to  allow  more  variability  in 
the  way  they  can  be  used.  Neill  cylinder  use  may  be  more  easily  reproduced,  which  is 
particularly  important  if  many  field  technicians  are  involved  in  the  collections.  The  time 
taken  to  obtain  five  Neill  samples  is  slightly  more  than  the  time  required  to  collect  a 
single  kick  sample,  but  not  significantly  so. 

Sample  processing  time  and  cost  are  also  a considération.  Even  if  the  cost  for  processing 
Neill  and  Kick  samples  are  similar,  processing  Neill  samples  would  be  costlier  because 
collection  of  replicates  is  the  norm.  Kick  samples  are  larger  than  individual  Neill  samples 
and  contain  more  débris  so  require  more  time  to  process  unless  severe  sub-sampling 
measures  are  taken.  Although  the  total  number  of  organisms  in  five  replicate  Neills  can 
be  fairly  close  to  the  total  found  in  a 210  pm  kick  sample,  more  of  the  material  tends  to 
get  sorted.  In  the  current  study,  even  the  coarse  2.0  mm  fraction  of  the  210  and  400  kicks 
from  Stier’s  Ranch  had  to  be  quartered  to  make  sorting  and  identifying  manageable.  For 
the  1 .0  mm  fraction,  most  of  the  Neill  replicates  were  sorted  in  their  entirety,  while  on 
average  half  of  the  fraction  was  sorted  for  the  210  kicks  and  a third  of  the  fraction  was 
sorted  for  the  400  kicks.  For  the  finest  fraction  (0.2  mm),  a little  more  than  a tenth  of  the 
Neill  material  and  a quarter  of  the  400  kick  sample  was  used  on  average.  For  ail  the  210 
kicks  except  from  Cochrane,  only  one  fortieth  of  the  fine  fraction  was  evaluated.  The 
smaller  the  sample  fraction  sorted  the  larger  the  uncertainty  surrounding  the  estimated 
abundance  and  richness  of  that  sample. 

3. 1.5  Conclusions 

Although  there  is  a higher  cost  associated  with  processing  Neill  cylinder  samples,  this 
method  appears  to  be  the  most  appropriate  for  continuing  to  monitor  longitudinal  and 
temporal  trends  in  the  benthos  of  Alberta’ s major  ri  vers.  The  method  is  repeatable  and 
relatively  simple.  The  replicates  collected  at  each  site  provide  smaller,  more  manageable 
samples  that  provide  greater  statistical  power  than  a single  kick  sample.  Considered 
together,  the  Neill  replicates  hâve  the  ability  to  reflect  greater  taxonomie  diversity  than 
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either  the  210  or  400  kick  samples.  The  quantitative  nature  of  the  Neill  cylinder  allows 
population  density  estimâtes  that  are  not  possible  with  kick  samples.  Furthermore,  given 
that  Neill  cylinder  samples  hâve  historically  been  used  in  Alberta  rivers  by  AENV,  the 
results  will  be  directly  comparable  with  20  years  of  previous  data. 

Many  of  the  site-to-site  différences  that  were  observed  in  the  Neill  samples  were 
reflected  in  the  kick  samples  as  well.  Kick  samples  may  very  well  be  appropriate  in  other 
contexts,  such  as  in  a régional  survey  of  smaller  streams  where  many  locations  are  being 
examined.  In  a case  where  kick  nets  are  preferred,  the  choice  of  mesh  size  will  be 
important.  Although  a smaller  mesh  allows  the  capture  of  more  organisms,  the  potential 
for  very  large  samples  that  may  be  biased  towards  smaller  organisms  must  be  considered. 

3.2  Assessment  Of  Longitudinal  And  Temporal  Trends  In  Benthic 
Invertebrate  Data 

3.2. 1 Longitudinal  Trends  in  Erosional  Habitat  (2006)  - Neill  Cylinder 
Samples 

3. 2. 1.1  Abundance , Richness , and  Diversitv 

Différences  in  the  numbers  and  types  of  invertebrates  found  at  each  of  the  five  Bow 
River  sites  show  longitudinal  patterns  consistent  with  natural  and  anthropogenic  changes 
along  the  length  of  the  river.  The  lowest  abundance  of  organisms  (total  number  of 
organisms)  was  found  at  Cochrane,  while  the  highest  was  found  at  Stier’s  Ranch  (Figure 
3.1).  Abundance  at  Cochrane  was  significantly  lower  than  at  ail  other  sites;  the  only 
other  significant  différence  was  between  Stier’s  and  Cluny  (Appendix  7).  The  Cochrane 
site  is  upstream  of  most  major  point-sources  of  nutrient  enrichment  on  the  Bow  River 
(except  for  some  municipal  and  industrial  discharge  from  the  Bow  Valley  Corridor).  In 
contrast,  Stier’s  Ranch  is  directly  below  the  City  of  Calgary  and  its  two  wastewater 
treatment  facilities,  which  provide  nutrients  to  maintain  substantial  populations  of  algae, 
macrophytes,  and  invertebrates. 

The  number  of  unique  taxa  (richness)  provides  information  about  the  “health”  or  degree 
of  perturbation  at  a site,  with  fewer  taxa  suggesting  a more  impacted  site.  In  the  Bow 
River,  the  lowest  richness  occurred  at  Cochrane  (Figure  3.2),  where  the  number  of  taxa 
was  significantly  lower  than  at  ail  other  sites  (Appendix  7).  Given  Cochrane’ s foothills 
location  and  heterogeneous  substrate,  the  site  would  be  expected  to  support  a larger 
variety  of  taxa.  However,  as  suggested  by  Anderson  (1991),  manipulation  of  flow  as  a 
resuit  of  the  operation  of  the  Ghost  River  hydroelectric  dam  upstream  of  the  site  causes 
significant  diurnal  disturbance  to  the  benthic  fauna  at  this  site.  Benthic  communities  in 
this  reach  of  the  river  may  be  submerged  and  exposed  within  a 24  hour  period,  depending 
on  their  location  in  the  channel.  This  situation  would  not  be  tolerated  by  many  sensitive 
organisms  and  would  affect  the  distribution  of  other  organisms.  Collecting  a 
représentative  sample  at  such  a site  is  difficult,  as  samplers  may  not  be  aware  that  they 
are  accessing  an  area  that  has  been  recently  exposed  to  air  (Anderson  1991). 

Furthermore,  during  periods  of  high  flows,  flow  velocity  and  steep  slopes  make  it 
difficult  to  penetrate  deep  enough  in  the  river.  Didymosphenia  geminata,  a colonial 
diatom,  was  encountered  regularly  when  collecting  invertebrates  in  the  1980’ s and 
1990’ s;  it  was  présent  in  2006  although  perhaps  not  to  the  same  degree.  The  diatom  can 
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cover  rocks  with  cotton  wad-like  material  and  contributes  further  to  a rather  inhospitable 
habitat  for  invertebrates.  The  distribution  of  Didymosphenia  geminata  in  Alberta  Rivers 
is  believed  to  be  linked  to  flow  régulation  (Kirkwood  et  al.  2007). 

For  richness,  as  for  ail  metrics  involving  counts  of  taxa,  the  Neill  cylinder  data  can  be 
summarized  as  either  the  mean  number  of  taxa  per  site  (based  on  5 replicates)  or  as  the 
total  number  of  unique  taxa  at  a site,  the  latter  better  reflecting  the  site’ s true  complément 
of  taxa,  including  rare  organisms.  In  this  case,  while  the  highest  mean  number  of  taxa 
occurred  at  Stier’s  Ranch,  the  highest  total  number  of  unique  taxa  occurred  at  Ronalane 
(Figure  3.2,  Appendix  7).  Note  also  most  of  the  analyses  for  the  2006  data  are  based  on  a 
“truncated”  data  set,  reflecting  taxonomie  levels  and  names  used  in  the  historical  data 
sets.  Using  the  original  taxon  list  for  the  2006  data,  total  richness  is  higher  at  each  site  by 
between  12  and  23  taxa  (Figure  3.3).  The  overall  pattern  is  similar  to  that  based  on 
higher  taxonomie  levels;  however,  the  highest  mean  and  total  numbers  of  taxa  occur  at 
Carseland,  mainly  due  to  a large  variety  of  chironomids  at  this  site. 

Because  they  summarize  the  relative  abundance  of  various  taxa,  it  is  expected  that 
diversity  indices  provide  more  succinct  information  about  benthic  communities  than 
abundance  or  richness  alone.  Two  such  common  indices,  Simpson’ s (JD)  and  Shannon’s 
(H)  were  calculated  for  the  2006  data.  The  Simpson  Index  showed  the  lowest  diversity  at 
Cochrane  and  Ronalane,  with  Cochrane  differing  significantly  only  from  Stier’s  Ranch, 
the  most  diverse  site  (Figure  3.4,  Appendix  7).  This  index  de-emphasizes  rare  taxa,  while 
highlighting  common  taxa  and  evenness  or  equability  among  taxa  (Mandaville  2002; 
Béais,  Gross,  and  Harrell  1999).  The  Shannon  Index  also  considers  abundance  and 
evenness  of  taxa,  but  is  more  affected  by  rare  taxa  (Mandaville  2002;  Béais,  Gross,  and 
Harrell  2000).  This  index  again  showed  Stier’s  Ranch  to  hâve  the  highest  diversity,  and 
Cochrane  the  lowest,  although  the  results  did  not  differ  significantly  (Figure  3.5, 
Appendix  7). 

3.2. 1.2  Counts  and  Provortions  of  Taxa 

The  previous  measures  suggest  that  Cochrane  expériences  some  disturbance  which  limits 
benthic  invertebrate  abundance  and  richness,  while  enrichment  downstream  of  Calgary 
supports  a larger  and  more  diverse  population.  More  spécifie  information  about  the  taxa 
at  each  site  is  necessary  to  describe  further  différences  in  the  benthic  communities  along 
the  length  of  the  Bow  River.  The  presence  of  a variety  of  Ephemeroptera,  Plecoptera, 
and  Trichoptera  (EPT)  taxa  is  generally  indicative  of  désirable  lotie  conditions.  The 
number  of  EPT  généra  was  significantly  lower  at  Cochrane  than  at  the  other  sites,  and 
was  also  significantly  lower  at  Ronalane,  near  the  mouth  of  the  Bow  River,  than  at  the 
previous  three  sites  (Figure  3.6,  Appendix  7).  Stoneflies  (Plecoptera)  were  completely 
absent  from  both  Cochrane  and  Ronalane.  While  flow  fluctuations  likely  limit  the 
number  of  EPT  taxa  at  Cochrane,  it  is  probable  that  many  of  these  taxa  are  excluded  from 
Ronalane  by  the  slower  flow  and  more  homogeneous  substrate  typical  of  a depositional 
(as  opposed  to  erosional)  environment.  Lower  dissolved  oxygen  levels  may  limit 
sensitive  taxa  at  Ronalane  as  well. 

In  erosional  habitats,  the  balance  between  the  abundance  of  EPT  organisms  and  the 
generally  more  pollution-tolérant  Chironomidae  should  be  relatively  even  in  a 
community  in  good  biotic  condition  (Mackie  2004).  The  EPT  to  Chironomidae  ratio 
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tended  to  be  highest  at  Stier’s  Ranch  (mean  0.59),  although  replicates  varied  widely, 
while  the  lowest  was  at  Ronalane  (mean  0.02)  (Figure  3.7).  The  mean  ratio  at  Stier’s 
Ranch  was  significantly  higher  than  that  at  the  Cochrane,  Carseland,  or  Ronalane  sites 
(Appendix  7). 

Chironomids,  nematodes  and  oligochaetes,  as  major  groups,  are  often  described  as  being 
tolérant  to  adverse  conditions,  especially  nutrient  enrichment  and  low  DO.  However, 
many  taxa  belong  to  these  groups  and  the  perceived  tolérance  may  be  attributable  to  only 
a few  taxa.  Unfortunately  this  is  seldom  recognized  because  species  or  genus 
identifications  are  rarely  performed  on  nematodes  and  oligochaetes.  The  highest 
proportions  of  chironomids  were  in  samples  taken  from  Carseland,  while  the  lowest 
proportions  were  found  at  Cochrane  (Figure  3.8).  Cochrane,  however,  had  the  highest 
proportion  of  Nematoda  (Figure  3.9).  Oligochaetes  were  présent  in  the  greatest 
proportion  at  Ronalane  folio wed  by  Cochrane  (Figure  3.10). 

The  relative  contribution  to  abundance  by  different  taxa  is  another  measure  of 
environmental  condition.  In  a less  complex,  potentially  more  impacted  site,  the  percent 
contribution  of  the  dominant  group  would  be  expected  to  be  higher  than  at  a less 
impacted  site.  In  terms  of  insects,  members  of  the  order  Diptera  were  most  abundant  at 
ail  5 sites.  The  percent  dominance  of  this  order  was  higher  at  Ronalane  (mean  98%)  than 
at  ail  other  sites,  and  significantly  so  for  ail  but  Carseland  (Figure  3.1 1,  Appendix  7). 

3.2. 1.3  Functional  Feeding  Groups 

Aquatic  invertebrates  can  be  classified  into  functional  feeding  groups  by  the  type  and  size 
of  food  particles  they  eat  and  by  the  method  that  food  is  obtained  (Merritt  and  Cummins 
1996).  The  feeding  group  concept  is  useful  for  describing  the  ecological  rôle  of  aquatic 
organisms.  The  numbers  of  invertebrates  with  various  feeding  strategies  can  then  be  used 
to  make  inferences  about  river  conditions. 

Filtering  collectors,  or  suspension  feeders,  feed  on  fine  particulate  organic  matter. 
Ronalane  had  the  highest  number  of  filter-feeding  taxa  (Figure  3.12,  Appendix  7). 
However,  the  proportion  of  filters  was  significantly  higher  at  Stier’s  Ranch  than  at  ail 
other  sites  (Figure  3.13,  Appendix  7),  likely  as  a resuit  of  organic  enrichment  from  the 
City  of  Calgary.  In  contrast,  gathering  collectors,  which  consume  coarser  particulate 
organic  matter,  occurred  in  the  highest  proportion  at  Ronalane  (Figure  3.14,  Appendix  7). 
Both  Cochrane  and  Stier’s  Ranch  had  much  lower  proportions  of  gatherers  than  did  the 
downstream  sites. 

Shredder  organisms,  which  “shred”  coarse  organic  matter  such  as  leaf  litter,  were  found 
in  the  highest  proportion  at  Cochrane  (Figure  3.15).  Shredder  proportions  were  relatively 
low  at  ail  sites,  however,  and  the  différences  between  the  sites  were  not  significant 
(Appendix  7).  Scraper  organisms,  which  feed  directly  on  periphyton,  were  found  in 
much  lower  proportions  at  Cochrane  and  Ronalane  than  at  the  other  sites  (Figure  3.16, 
Appendix  7),  reflecting  lower  amounts  of  epilithic  algae  at  these  sites. 

Predators,  which  feed  on  other  invertebrates,  were  found  in  the  highest  proportion  at 
Stier’s  Ranch  (mean  0.10),  the  site  with  the  highest  invertebrate  abundance.  Predator 
proportions  were  the  lowest  at  Ronalane  (mean  0.03)  (Figure  3.17,  Appendix  7). 
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3.2. 1.4  Habits  and  Habitat  Preference 


Différences  in  the  behaviour  or  habit  of  benthic  invertebrates,  defined  as  the  way  they 
move  or  maintain  position  in  flowing  water,  were  examined  for  the  Bow  River  sites.  The 
proportion  of  swimmers  was  significantly  higher  at  Cochrane  than  anywhere  else  (Figure 
3.18,  Appendix  7).  Swimming  taxa  hâve  the  greatest  motility,  and  thus  would  be  best 
able  to  deal  with  the  dynamic  flow  and  water  level  fluctuations  at  this  site. 

Clinger  organisms  hold  their  position  on  the  bottom  substrate  in  flowing  water.  Greater 
numbers  and  proportions  of  clinger  taxa  are  expected  under  conditions  subject  to  less 
impact  from  physical  disturbance  or  pollution.  Cochrane  and  Ronalane  had  significantly 
fewer  clinger  taxa  than  the  other  sites  (Figure  3.19,  Appendix  7).  In  terms  of  proportion, 
Stier’s  Ranch  had  significantly  more  clingers  and  Ronalane  significantly  less,  than  the 
other  sites  (Figure  3.20,  Appendix7).  Sprawler  taxa,  which  hold  on  to  both  the  substrate 
and  submerged  végétation,  showed  a similar  pattern,  with  the  highest  mean  proportion 
found  at  Stier’s  and  the  lowest  at  Ronalane  (Figure  3.21,  Appendix  7).  In  contrast,  the 
greatest  number  of  burrowing  taxa,  which  live  in  the  sédiment,  occurred  at  Ronalane 
(Figure  3.22).  The  proportion  of  burrowers  was  lowest  at  Stier’s  Ranch,  and  highest  at 
Ronalane,  where  on  average  95%  of  the  organisms  présent  were  burrowers  (Figure  3.23, 
Appendix  7). 

Given  the  pattern  shown  for  benthic  invertebrate  habits,  it  would  be  expected  that  habitat 
preference  (whether  erosional  or  depositional)  would  also  be  a useful  metric  for 
describing  longitudinal  patterns  in  the  Bow  River.  However,  based  on  the  taxonomie 
units  used  for  comparison  with  the  historical  data,  habitat  preference  was  misleading. 
There  were  no  interprétable  différences  among  sites  based  on  the  mean  proportion  of 
organisms  preferring  depositional  habitats;  Ronalane,  which  is  clearly  a slower-flowing 
depositional  site,  showed  the  lowest  proportion  of  depositional  organisms  (Figure  3.24, 
Appendix  7).  However,  when  taxonomie  detail  was  added  back  in,  Ronalane  stood  out 
from  the  other  sites  as  the  location  where  the  greatest  proportion  of  the  benthos  preferred 
depositional  habitats,  consistent  with  the  nature  of  this  site  (Figure  3.  25).  This 
différence  was  driven  by  the  spécifie  preferences  of  various  chironomid  généra. 

3.2. 1.5  Tolérance 

The  tolérance  or  sensitivity  of  many  benthic  invertebrates  to  pollution  and  other 
perturbation  is  fairly  well  known.  For  this  report,  pollution  tolérance  and  pollution 
sensitivity  grades  were  assigned  to  taxa  based  on  tables  in  Mandaville  (2002),  which 
were  drawn  from  a number  of  original  sources.  Tolérance  values,  whether  at  the 
species/genus  level  or  the  family  level,  range  from  0 (very  intolérant)  to  10  (very 
tolérant)  and  are  based  primarily  on  response  to  organic  pollution.  These  values  can  be 
combined  into  a biotic  index  (BI)  or  a family  biotic  index  (FBI)  to  give  a range  of  values 
between  0 and  10,  where  0 represents  excellent  water  quality  and  10  represents  very  poor 
water  quality.  Pollution  sensitivity  grades  are  assigned  on  a family  basis  or  higher. 

These  values  range  from  1 to  10,  where  1 is  least  sensitive  and  10  is  most  sensitive.  The 
sensitivity  grades  can  be  examined  as  the  cumulative  score  of  ail  taxa  (referred  to  as  the 
Biological  Monitoring  Working  Party  number,  BMWP)  or  can  be  divided  by  the  number 
of  taxa  to  yield  an  average  score  per  taxon  (ASPT). 
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The  number  of  pollution  intolérant  taxa  (score  of  0,  1 , or  2)  was  lowest  at  Ronalane, 
followed  by  Cochrane  (Figure  3.26,  Appendix  7).  There  were  no  significant  différences 
between  the  other  three  sites,  although  Carseland  had  the  highest  mean  number  of 
intolérant  taxa,  and  was  tied  with  Cluny  for  the  highest  total  number  of  intolérant  taxa. 

The  mean  Family  Biotic  Index  (FBI)  scores  for  the  5 Bow  River  sites  ranged  from  7.8 
(Cochrane)  to  5.8  (Stier’s  Ranch).  Scores  at  Cochrane  and  Ronalane  were  significantly 
higher  than  at  the  other  three  sites  indicating  that  these  two  sites  sustained  the  most 
négative  impact  (Figure  3.27,  Appendix  7).  Although  schemes  exist  for  translating  BI 
and  FBI  scores  into  water  quality  ratings  ( e.g . Hilsenhoff  1987,  1988),  the  values 
calculated  for  the  Bow  sites  should  be  interpreted  as  relative  rankings  only.  Attributing 
the  scores  to  spécifie  water  quality  ratings  can  only  be  done  after  extensive  calibration  for 
the  région  in  which  they  are  to  be  used. 

Mean  cumulative  sensitivity  (BMWP)  scores  were  significantly  lower  at  Cochrane  than 
ail  other  sites  except  for  Ronalane  (Figure  3.28,  Appendix  7).  The  general  paucity  of 
fauna  at  the  Cochrane  site  had  some  bearing  on  this  resuit.  Once  the  effect  of  the  number 
of  taxa  was  eliminated  by  calculating  the  average  sensitivity  score  per  taxon,  the  gap 
between  Cochrane  and  the  other  sites  was  narrowed,  and  the  significance  of  the 
différence  between  Cochrane  and  Stier’s  Ranch  was  eliminated.  Both  Cochrane  and 
Ronalane  had  significantly  lower  ASPT  than  Cluny  and  Carseland  (Figure  3.29, 

Appendix  7).  These  results  suggest  that  although  nutrient  enrichment  at  Stier’s  Ranch 
resulted  in  a rich  and  abundant  benthic  community,  there  may  be  a slight  tendency  for 
those  organisms  to  be  less  sensitive  to  organic  pollution  than  those  in  adjacent 
downstream  communities,  farther  from  Calgary’s  direct  influence. 


3.2.2  Longitudinal  t rends  in  Depositional  Habitat  (2006)  - 210  pm 
Kick  Net  Samples 

Most  of  the  benthic  invertebrate  monitoring  programs  by  Alberta  Environment  hâve 
focused  on  erosional  habitat  and  the  benthic  data  collected  from  depositional  areas  at 
Bow  River  LTRN  sites  are  among  the  first  for  such  habitat. 

In  many  ways  the  longitudinal  patterns  exhibited  by  benthic  invertebrate  metrics  in 
depositional  areas  mirror  those  observed  in  erosional  areas.  For  example,  the  number  of 
invertebrates  collected  downstream  of  Calgary  was  considerably  higher  than  upstream 
(Figure  4.1),  indicating  that  responses  to  enrichment  are  also  measurable  in  depositional 
areas. 

The  number  of  taxonomie  groups  collected  in  depositional  areas  tended  to  be  lower  than 
in  erosional  areas  sampled  with  the  same  mesh  size  (Figure  4.2).  This  is  particularly 
noticeable  for  some  taxa  such  as  mayflies,  stoneflies,  and  caddisflies  which  were 
considerably  less  diverse  and  less  abundant  in  depositional  areas  (Figures  4.3  and  4.4). 
As  expected,  the  proportion  of  taxa  which  are  known  to  prefer  depositional  habitats  was 
greater  in  depositional  samples  than  in  erosional  samples  (Figure  5)  and  taxa  which 
burrow  in  sédiments  were  more  prévalent  in  this  habitat  (Figure  6).  Swimmers,  which 
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are  quite  mobile  (over  the  sédiment  and  in  the  water  column),  also  tended  to  be  relatively 
more  abundant  in  depositional  areas. 

Unlike  erosional  habitats  where  diversity  indices  change  relatively  little  along  the  length 
of  the  Bow  River,  there  was  a fairly  sharp  longitudinal  décliné  in  index  values  in 
depositional  habitats  downstream  of  Calgary  and  recovery  was  only  apparent  at  the  most 
downstream  site  (Figure  4.5  and  4.6).  By  nature,  depositional  areas  hâve  slower  flow 
velocities.  Dissolved  oxygen  fluctuations  may  be  more  pronounced  than  in  faster  flowing 
waters,  particularly  in  areas  that  are  very  rich  in  nutrients,  and/or  hâve  high  algal  and 
macrophyte  biomass.  Subséquent  décomposition  of  primary  producers  may  lead  to 
excessive  oxygen  consumption  from  the  water  column.  The  décliné  in  diversity  indices 
may  be  an  indication  of  more  stressful  conditions  in  depositional  than  erosional  areas.  If 
this  can  be  verified,  it  may  be  judicious  to  continue  monitoring  such  depositional  areas  to 
assess  responses  to  nutrient  management  in  the  basin. 


3.2.3  Historical  Longitudinal  and  Temporal  Trends  for  Invertebrates 
in  Erosional  Habitat 

Many  of  the  patterns  described  from  the  2006  data,  including  evidence  of  perturbation  at 
Cochrane,  organic  enrichment  at  Carseland,  and  depositional  conditions  at  Ronalane  are 
consistent  with  the  longitudinal  différences  in  the  Bow  River  benthic  community 
discussed  in  Anderson  (1991).  This  historical  work  was  based  on  data  collected  between 
1983  and  1987  at  Cochrane,  Carseland  and  Ronalane.  For  the  current  report,  these  data 
were  combined  with  occasional  data  from  Stier’s  Ranch  and  Cluny  and  data  collected  in 
1989,  1993,  and  2006  in  order  to  examine  overall  différences  among  sites  and  over  time 
using  a number  of  different  techniques,  including  multivariate  analyses. 

3.2.3. 1 Steywise  Discriminant  Analysis  - Ail  sites  and  Years 

Stepwise  discriminant  analysis  was  performed  using  counts  of  the  top  25  taxa  (out  of  a 
total  of  64).  Four  canonical  discriminant  functions  were  used  in  the  analysis,  the  first 
three  of  which  accounted  for  about  96%  of  the  variance  among  samples.  The  corrélations 
for  these  functions  are  shown  in  Appendix  8.  Of  the  25  taxa  used  in  the  discriminant 
analysis,  14  were  associated  with  the  first  3 functions. 

Samples  from  Cochrane  were  clearly  separated  from  the  rest  of  the  Bow  River  sites  along 
function  1 , based  primarily  on  fewer  Chironomini,  Tubificidae,  and  Simulidae  at  this  site 
than  others  (Figure  7).  Ail  sites  showed  more  séparation  along  function  2,  which  was 
positively  correlated  with  the  mayflies  Stenonema  and  Tricorythodes,  and  negatively  with 
Hydra  and  Ephemerellidae.  Cochrane  and  Stier’s  Ranch  held  a neutral  position  along 
this  axis,  while  Carseland  had  a négative  score  and  both  Ronalane  and  Cluny  had  positive 
scores  for  this  function  (Figure  7).  Discriminant  function  3 was  positively  correlated 
with  the  caddis  fly  family  Helicopsychidae  and  to  a lesser  extent  with  the  chironomid 
tribe  Tanytarsini.  Carseland,  Cluny,  and  Ronalane  scored  higher  on  this  function  and 
Cochrane  occasionally  scored  lower,  while  Stier’s  Ranch  samples  scored  negatively 
(Figure  7).  Considering  ail  three  functions,  Ronalane  and  Cluny  were  the  most  similar, 
showing  the  most  overlap,  and  the  other  three  sites  each  formed  distinct  separate  groups 
(Figure  7). 
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The  fairly  good  séparation  of  the  five  sites  based  on  this  analysis  suggests  that  many  of 
the  basic  longitudinal  différences  in  the  zoobenthic  community  of  the  Bow  River 
remained  consistent  over  the  23  year  period  of  study.  However,  to  better  detect 
significant  trends  over  time  within  sites,  analysis  of  variance  (ANOVA)  was  performed 
using  the  discriminant  function  scores.  The  hypothesis  was  that  no  change  had  occurred 
in  the  taxa  influencing  the  function  in  question,  while  the  alternate  hypothesis  was  that 
significant  change  had  occurred.  For  DF-1,  only  Stier’s  Ranch  showed  a significant 
increasing  change  over  time  (Figure  8).  The  taxa  most  strongly  associated  with  this 
function  tend  to  be  pollution  tolérant,  suggesting  a trend  towards  more  tolérant  organisms 
at  this  site.  For  DF-2,  scores  at  Ronalane  decreased  significantly  over  time.  Taxa 
associated  with  this  function  were  less  tolérant  to  pollution  than  those  for  DF-1, 
suggesting  a trend  towards  fewer  sensitive  organisms  at  this  site.  For  DF-3,  Carseland’s 
scores  significantly  decreased  over  time,  while  Cluny’s  scores  increased  significantly. 
DF-3  was  most  strongly  associated  with  the  presence  of  helicopsychid  caddisflies  which 
prefer  erosional  habitats.  Overall,  these  analyses  suggest  some  change  over  time  at  ail 
sites  except  for  Cochrane.  At  ail  sites,  it  appears  that  conditions  may  hâve  been  more 
favourable  to  pollution-sensitive  organisms  in  the  middle  of  the  period  studied  than  in  the 
early  eighties  or  at  présent,  with  some  minor  exceptions. 

3.2. 3. 2 Comyarison:  1985 , 1993  and  2006 

There  are  only  three  years  with  data  at  ail  five  sites:  1985,  1993,  and  2006.  These  three 
years  represent  fairly  well-spaced  time  periods:  post-phosphorus  removal  (1980s),  post- 
nitrogen  removal  (1990s),  and  présent  (2000s). 

Réductions  in  nutrient  loading  would  be  expected  to  resuit  in  less  primary  production  and 
improved  dissolved  oxygen  conditions.  Depending  on  the  degree  of  these  changes, 
benthic  invertebrate  communities  would  be  expected  to  show  increased  taxonomie 
richness,  particularly  for  ‘sensitive’  taxa.  If  nutrient  réductions  were  sufficient  to  resuit 
in  a décliné  of  invertebrate  food  resources,  a décliné  in  overall  population  densities  might 
also  be  seen. 

The  flow  history  that  précédés  invertebrate  sampling  can  hâve  a determining  influence  on 
benthic  invertebrate  distribution  and,  in  this  context,  an  évaluation  of  flows  has  been 
included  in  Appendix  9.  Flow  conditions  were  quite  different  in  the  three  years.  Flows 
preceding  the  sampling  in  1985  were  generally  low,  but  a sharp  peak  occurred  a month 
earlier  in  September.  Flows  for  the  2006  sampling  were  higher,  but  had  remained  fairly 
stable  for  much  of  the  summer  and  fall;  very  high  flows  occurred  in  June  2006.  The 
September  1993  samples  represent  the  highest  monthly  and  daily  flows  of  ail  years  and 
numerous  high  flow  events  occurred  during  the  summer.  The  higher  flows  experienced 
in  1993,  in  many  ways  could  hâve  favoured  ‘sensitive’  taxa. 

A graphical  comparison  of  benthic  invertebrate  metrics  for  1985,  1993  and  2006  follows. 
This  comparison  is  supported  by  a sériés  of  ANOVAs,  examples  of  which  are  provided  in 
Appendix  10). 

3.2. 3. 3 Counts  and  Indices 
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In  spite  of  some  site-to-site  différences,  overall  abundance  (measured  as  the  mean  of  the 
five  sites)  was  similar  for  1993  and  2006,  while  it  was  much  lower  in  1985  (Figure  9-1). 

It  is  possible  that  high  fiows  prior  to  sampling  resulted  in  lower  benthos  numbers  in  that 
year.  However,  in  terms  of  richness  (shown  for  number  of  families),  1985  and  2006  were 
similar,  while  the  number  of  taxa  was  generally  higher  in  1993,  which  could  be  an 
indication  of  more  favourable  conditions  (Figure  9-2).  Both  abundance  and  richness 
were  significantly  lower  at  Cochrane  in  2006  than  in  the  previous  years.  At  Stier’s 
Ranch,  abundance  was  much  lower  in  1985  than  in  2006,  while  the  number  of  families 
found  in  these  two  years  was  not  significantly  different. 

Although  the  overall  Simpson  diversity  index  was  highest  in  1993,  the  only  significant 
différence  between  years  was  at  Cluny,  which  had  a much  less  diverse  benthic 
community  in  1985  than  in  the  two  other  years  (Figure  9-3).  The  overall  Shannon 
diversity  index  was  highest  in  1993  (Figures  9-4).  The  2006  index  value  at  Cochrane  was 
significantly  lower  than  it  had  been  in  1993,  and  at  Ronalane  the  2006  index  was 
significantly  lower  than  in  both  other  years.  The  Shannon  index  was  most  variable  at 
Cluny;  ail  three  years  were  different,  with  1985  having  the  lowest  value  and  2006  the 
highest. 

Greater  community  richness  throughout  the  Bow  River  in  1993  was  reflected  in  overall 
higher  numbers  of  taxa  (généra  or  families)  in  various  categories,  whether  based  on 
feeding  groups,  behaviour,  or  pollution  sensitivity.  However,  the  patterns  at  individual 
sites  varied.  For  Ephemeroptera,  Plecoptera,  and  Trichoptera  (EPT)  généra,  numbers 
were  significantly  higher  in  1993  at  ail  sites  than  in  2006,  and  higher  at  Stier’s  Ranch  and 
Cluny  in  1993  than  in  1985  (Figure  9-5).  At  Cochrane  and  Ronalane,  the  number  of 
généra  of  these  three  “sensitive”  groups  was  also  larger  in  1985  than  in  2006.  In  terms  of 
families  within  the  orders  Ephemeroptera,  Trichoptera,  and  Odonata  (ETO,  another 
metric  for  grouping  sensitive  organisms),  there  was  no  différence  among  years  for  Stier’s 
Ranch  or  for  Carseland,  while  counts  in  1985  and  1993  were  significantly  higher  than 
those  in  2006  for  both  Cochrane  and  Ronalane  (Figure  9-6).  Numbers  for  Cluny  were 
significantly  higher  in  1993  than  in  both  other  years. 

Numbers  of  various  taxa  at  Cochrane  often  showed  the  largest  différences,  particularly 
between  1985  and  2006.  For  example,  the  largest  number  of  predator  taxa  for  ail  sites 
and  years  was  found  at  Cochrane  in  1985,  while  the  lowest  number  occurred  at  the  same 
site  in  2006  (Figure  9-7).  At  Carseland  and  Cluny,  the  highest  number  of  predator  taxa 
occurred  in  1993,  contributing  to  high  overall  numbers  at  this  site,  while  at  Stier’s  Ranch 
and  Ronalane,  there  were  no  significant  différences  among  years.  The  pattern  for  clinger 
taxa  was  similar;  the  highest  overall  numbers  occurred  in  1993,  while  the  highest  and 
lowest  number  of  clinger  taxa  occurred  at  Cochrane  in  1985  and  2006,  respectively 
(Figure  9-8).  The  Cochrane  site  was  not  affected  by  the  same  September  peak  flow  that 
the  downstream  sites  experienced  before  sampling  in  1985  (Appendix  9),  which  may 
hâve  displaced  taxa  at  these  sites.  Lower  flows  in  1985  may  also  hâve  allowed  samplers 
to  access  an  area  of  the  river  bed  at  Cochrane  that  was  not  as  much  disturbed  by 
discharge  fluctuations  resulting  from  operation  of  the  Ghost  Dam,  and  had  a more 
représentative  fauna  than  the  2006  samples. 

The  overall  number  of  pollution  intolérant  taxa  was  almost  twice  as  high  in  1993  as  in  the 
other  two  years.  In  1985,  pollution  intolérant  taxa  were  more  common  at  Cochrane  than 
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at  any  other  site,  but  in  2006,  the  number  of  intolérant  taxa  at  this  site  was  second  only  to 
Ronalane’s  count  that  year  as  the  lowest  for  ail  years  (Figure  9-9).  The  only  site  that  did 
not  differ  significantly  among  years  was  Carseland.  This  was  the  case  also  for  the 
average  pollution  sensitivity  score  per  taxon,  which  showed  similar  year-to-year  and  site- 
to-site  patterns  (Figure  9-10). 

The  Biotic  Index  (based  on  data  at  the  genus  level  or  higher)  produces  scores  between  1 
and  10  where  higher  scores  correspond  to  poorer  water  quality.  The  highest  biotic  index 
scores  occurred  at  Carseland  and  Cluny  in  1985;  these  were  significantly  higher  (poorer) 
than  in  1993  or  2006.  At  Cochrane,  the  index  score  in  2006  was  significantly  higher  than 
in  the  previous  years.  At  Stier’s  and  Ronalane  différences  among  years  were  not 
statistically  significant  (Figure  9-11).  As  for  most  other  metrics,  the  overall  biotic  index 
results  were  lowest  (best)  in  1993.  High  flows  throughout  the  summer  may  hâve 
constantly  refreshed  the  substrate  and  kept  oxygen  levels  high  in  that  year  (Appendix  9). 
In  1993  both  phosphorus  and  nitrogen  removal  were  fully  implemented,  but  the 
population  of  Calgary  was  less  than  three  quarters  of  what  it  was  in  2006  (City  of  Calgary 
2007).  The  alum  treatment  and  resulting  effluent  phosphorus  réduction  may  hâve  helped 
to  offset  the  loading  increase  due  to  population  growth  (A.  Sosiak,  pers.  comm.). 

3.2. 3.4  Proportions  and  Ratios 

The  proportion  of  Oligochaeta  in  the  samples  differed  little  at  each  site  over  time,  except 
for  being  lower  in  2006  than  1993  at  Stier’s  Ranch,  and  lower  in  1985  than  1993  at 
Ronalane  (Figure  9-12).  The  proportion  of  Nematoda  differed  more  appreciably  among 
years,  with  the  overall  proportion  being  much  lower  in  1993  samples  than  in  other  years 
(Figure  9-13).  The  overall  proportion  of  Chironomidae  was  also  lower  in  1993  samples 
(Figure  9-14),  although  différences  at  Stier’s  Ranch  and  Carseland  were  not  statistically 
significant.  The  overall  proportion  of  EPT  organisms  and  the  corresponding  EPT-to- 
Chironomidae  ratio  showed  the  opposite  trend;  these  two  measures  were  much  higher  in 
1993  than  in  the  other  years  (Figures  9-15  and  9-16).  Stier’s  Ranch  had  slightly  higher 
proportions  of  EPT  compared  to  the  other  sites  in  ail  years,  except  for  Cluny  in  1993. 
Cluny  had  one  of  the  lowest  EPT-to-chironomid  ratios  in  1985  (mean  of  0.1 1),  while  in 
1993  the  mean  ratio  at  this  site  was  6.25  (Figure  9-16). 

The  percent  by  which  one  type  of  organism  dominâtes  in  a sample  provides  an  indication 
of  balance  in  the  community.  The  contribution  of  the  dominant  order  of  insects  at  each 
site  as  a percentage  of  total  insect  abundance  was  generally  higher  in  1985  and  2006.  In 
these  two  years,  maximum  contributions  at  some  sites  were  upwards  of  90%,  while  in 
1993,  the  maximum  percent  contribution  was  78%  at  Cluny  (Figure  9-17). 

The  ratio  of  relatively  tolérant  hydropsychid  caddisflies  to  total  Trichoptera  tends  to 
increase  in  response  to  organic  enrichment.  Hydropsychids  are  net  spinning  (filter 
feeding)  caddisflies  which  are  one  of  the  most  notable  invertebrates  at  Stier’s  Ranch 
because  of  their  abundance  and  size.  The  ratio  of  hydropsychids  to  caddisflies  at  this  site 
was  somewhat  lower  in  1993  than  in  the  two  other  years  (Figure  9-18).  At  Carseland  that 
ratio  was  considerably  higher  in  2006. 

The  longitudinal  différences  in  the  proportions  of  filtering  collectors  was  preserved 
between  years,  with  Stier’s  Ranch  having  the  greatest  proportion  of  filterers  in  each  year 
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(Figure  9-19).  Proportions  of  gathering  collectors  were  always  much  higher  than  those 
for  filtering  collectors;  this  metric  tended  to  be  highest  in  ail  years  at  Ronalane  (Figure  9- 
20). 

The  ratio  of  scrapers,  which  prefer  areas  with  a large  amount  of  algae,  to  filterers,  which 
feed  on  fine  particulate  organic  matter  and  are  sensitive  to  any  toxic  material  bound  to 
these  particles,  can  provide  information  about  aquatic  conditions.  This  ratio  changed 
noticeably  at  some  sites  between  years,  particularly  at  Ronalane  and  Cluny  (Figure  9-21), 
but  did  not  appear  to  follow  any  particular  pattern,  nor  was  there  any  association  between 
this  metric  and  algal  data. 

The  overall  proportion  of  clinger  taxa,  which  grasp  onto  rocks,  plants,  and  débris,  was 
greatest  in  1993.  In  each  year,  the  mean  proportion  of  the  invertebrate  community  made 
up  by  clingers  was  highest  at  Stier’s  Ranch,  although  in  1993  the  maximum  proportion 
was  found  at  Cochrane  (Figure  9-22).  In  contrast,  the  overall  proportion  of  animais  that 
burrow  into  the  substrate  was  lowest  in  1993,  perhaps  reflecting  the  fact  that  high  flows 
kept  the  river  bed  more  clear  of  sédiment  in  that  year.  The  largest  différence  between 
years  was  seen  at  Cluny  where,  in  1985,  burrowers  comprised  87%  of  the  organisms 
présent,  while  in  1993  only  14%  of  the  organisms  were  burrowers  (Figure  9-23).  This 
différence,  along  with  others  such  as  the  increase  in  EPT  organisms  at  this  site  between 
1985  and  1993  indicate  that  the  nature  of  this  site  changed  markedly  in  the  intervening 
years,  perhaps  due  to  displacement  of  sédiment  from  the  river  bed  and/or  decreased  algal 
biomass  as  a resuit  of  flushing  flows  (see  Appendix  9).  Another  possibility  is  that  there 
were  différences  in  the  exact  location  from  where  samples  were  collected. 

Changes  in  the  condition  of  a site  between  years  could  be  revealed  further  by  an 
examination  of  habitat  preference.  For  example,  it  appears  Cluny  may  hâve  changed 
from  a depositional  site  to  an  erosional  site  between  1985  and  1993.  However,  as 
indicated  in  the  previous  analysis  of  the  2006  data,  many  of  the  invertebrate 
identifications  for  the  historical  data  are  not  spécifie  enough  for  adéquate  assignment  of 
habitat  preference. 


3.3  Assessment  Of  Longitudinal  Trends  In  Epilithic  Algal  Data 
3.3. 1 Epilithic  Algae  - 2006 

3. 3. 1.1  Community  Analysis 

A total  of  104  epilithic  algae  taxa,  the  majority  of  which  were  diatoms,  were  identified 
from  the  five  Bow  River  sites  in  October  of  2006  (Appendix  3).  The  greatest  algal 
density  (units/cm2)  was  recorded  at  Carseland  and  was  primarily  due  to  diatoms  (Figure 
10-1).  Diatoms  made  up  the  greatest  portion  of  density  at  ail  sites.  Cyanobacteria  were 
the  second  largest  contributor  to  density  at  ail  sites  except  for  Cluny,  where 
Chlorophyceae  was  the  second  densest  group.  The  highest  density  of  cyanobacteria 
occurred  at  Stier’s  Ranch  (Figure  10-1). 

Diatoms  made  up  the  majority  of  the  biomass  (pg/cm2,  estimated  based  on  cell  volume) 
at  ail  sites;  biomass  calculated  by  this  method  was  much  higher  at  Carseland  than  at  the 
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other  sites  (Figure  10-2).  Stier’s  Ranch  and  Carseland  had  the  highest  biomasses  of 
cyanobacteria,  while  Cochrane  and  Ronalane  had  the  highest  biomasses  of 
Chlorophyceae.  Higher  density  and  biomass  of  cyanobacteria  is  consistent  with  organic 
enrichment,  while  Chlorophyceae  can  be  indicative  of  cleaner  water  sites.  At  Stier’s 
Ranch,  cyanobacteria  made  up  23%  of  the  total  biomass,  while  at  Carseland, 
cyanobacteria  contributed  only  3%  of  the  total  biomass  (Figure  10-3). 

The  fewest  taxa  were  found  at  Cochrane  (24),  while  Ronalane  had  the  greatest  taxa 
richness  (50)  (Figure  10-4).  In  terms  of  unique  généra,  richness  was  again  lowest  at 
Cochrane  (15),  while  Carseland  had  the  highest  number  (22)  (Figure  10-5).  Taxa  were 
used  to  calculate  the  Shannon  index  at  each  site,  which  followed  the  general  pattern  for 
taxon  richness  (Figure  10-6). 

Information  on  algal  community  composition  can  be  used  in  a manner  similar  to  the 
benthic  invertebrate  community  for  assessing  environmental  conditions.  Although  algal 
classes  broadly  suggest  whether  conditions  are  relatively  “clean”  or  subject  to  nutrient 
enrichment,  taxa  within  each  class  exhibit  a range  of  responses  to  nutrient  levels.  These 
properties  hâve  been  used  to  develop  a number  of  water  quality/nutrient  indices, 
including  some  based  on  the  taxonomie  composition  of  cyanobacteria  (Douterelo, 

Perona,  and  Mateo  2004)  and  diatoms  (Kelly  1998,  Jüttner  et  al.  2003,  Potapova  and 
Charles  2007).  The  dominance  of  diatoms  at  ail  the  Bow  River  sites  make  these  algae 
attractive  as  indicator  taxa.  Based  on  data  from  U. S.  rivers,  Potapova  and  Charles  (2007) 
compiled  a list  of  diatom  species,  indicating  whether  they  were  associated  with  low  or 
high  total  phosphorus  (<  10  pg/L,  > 100  pg/L)  and  with  low  or  high  total  nitrogen  (<  0.2 
pg/L,  >3.0  pg/L).  Since  these  concentration  ranges  are  applicable  to  the  Bow  River 
(e.g.,  Sosiak  2002),  the  broad  preference  ranges  were  used  at  a scoping  level  to  classify 
Bow  River  diatom  species  simply  as  either  ‘low’  or  ‘high’  TP  species,  and  ‘low’  or  ‘high’ 
TN  species.  Actual  indicator  values  provided  in  Potapova  and  Charles  (2007)  were  not 
used  in  this  analysis. 

By  number  of  taxa  and  by  biomass,  the  greatest  proportion  of  diatoms  preferring  both 
low  TP  and  TN  (unclassified  taxa  excluded)  were  found  at  Cochrane,  consistent  with  this 
site’ s location  upstream  from  major  sources  of  enrichment.  Carseland  had  the  next 
highest  proportion  of  diatom  taxa  and  biomass  preferring  lower  nutrients  (Figures  11-1  to 
1 1-4).  The  greatest  densities  of  low  nutrient-preferring  taxa,  however,  appeared  at  Cluny 
and  Ronalane.  At  Ronalane,  the  densities  of  nutrient-tolerant  and  intolérant  taxa  were 
equal,  while  at  ail  other  sites  the  tolérant  taxa  were  found  at  greater  densities  (Figures  1 1- 
5 and  11-6). 


Chlorophvll-a  and  Phaeophytin-a 

Chlorophyll-a  pro vides  a reliable  indicator  of  algal  biomass;  phaeophytin-a  is  a 
breakdown  product  of  chlorophyll-a.  Between  3 and  9 samples  for  chlorophyll-a  and 
phaeophytin-a  were  collected  at  each  site.  At  Stier’s  Ranch,  Carseland,  and  Cluny,  more 
than  one  person  collected  the  samples,  adding  to  the  range  of  variation  seen  at  these  sites 
and  illustrating  the  value  of  consistency  in  field  staff  and  training  for  spatial  and  temporal 
évaluations  (Figure  12-1).  In  each  case,  the  samples  hâve  been  pooled  for  analysis.  The 
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highest  mean  chlorophyll-Æ  (299  mg/m2)  was  found  at  Stier’s  Ranch,  while  the  lowest 
was  18.46  mg/m2  at  Cochrane  (Figure  12-2).  Mean  phaeophytin-Æ  was  highest  at 
Carseland,  indicating  a high  amount  of  turn-over  at  this  site,  and  suggesting  that  a large 
amount  of  dead  cells  contributed  to  the  high  density  and  biomass  noted  by  the  taxonomie 
analysis  at  this  site  (Figure  12-3).  The  chlorophyll-to-phaeophytin  ratio  was  highest  at 
Stier’s  Ranch,  indicating  rapid  growth  at  this  site.  This  ratio  was  lowest  at  Cochrane  and 
Ronalane,  suggesting  a slower  rate  of  growth  in  comparison  to  algal  décliné  (Figure  12- 
4). 


Corrélations  with  Invertebrate  Data 


Epilithic  algal  community  metrics,  including  density,  biomass,  taxa  richness,  and  the 
Shannon  index,  along  with  chlorophyll-rz  data,  were  compared  to  the  2006  invertebrate 
abundance,  taxa  richness,  and  Shannon  index  data  using  simple  corrélation.  The 
strongest  relationships  were  between  chlorophyll-a  and  invertebrate  abundance  (Figure 
13-1,  r = 0.87),  the  chlorophylle  to  phaeophytine  ratio  and  the  invertebrate  Shannon 
index  (Figure  13-2,  r = 0.85),  and  algal  richness  and  invertebrate  richness  (Figure  13-3,  r 
= 0.85).  These  relationships  are  not  unexpected  since  periphyton  is  a major  food  source 
for  many  invertebrate  taxa. 

3.3.2  Comparison  with  Historical  Data 

Limited  historical  epilithic  chlorophylle  data  are  available  for  most  of  the  Bow  River 
sites.  Data  from  2006  were  compared  with  data  collected  in  September  of  1993  (three 
replicates  per  site)  and  pooled  data  for  autumn  samples  from  1983  through  1986.  Data 
were  pooled  to  represent  conditions  in  the  mid-eighties  (after  phosphorus  removal  but 
before  nitrogen  réduction)  because  only  one  sample  was  collected  for  each  sampling 
month  during  this  period.  No  fall  data  are  available  for  Cochrane  from  this  period. 

The  highest  mean  and  maximum  chlorophylle  levels  occurred  at  Carseland  and  Cluny  in 
the  1980s.  Algal  biomass  at  these  two  sites  was  lower  in  1993  and  2006.  In  both  1993 
and  2006,  the  largest  mean  Chie  values  were  found  at  Stier’s  Ranch,  followed  by 
Carseland  (Figure  14).  These  results  are  consistent  with  those  of  Sosiak  (2002),  who 
found  that  periphytic  biomass  did  not  change  directly  downstream  from  Calgary  in 
response  to  enhanced  phosphorus  removal,  but  declined  at  sites  further  downstream 
where  total  dissolved  phosphorus  levels  were  lower. 

The  chlorophyll  results  relate  well  to  the  invertebrate  findings  for  these  three  periods.  In 
particular,  the  lower  amounts  of  algae  in  1993,  perhaps  in  part  a resuit  of  higher  flows 
throughout  the  summer,  are  consistent  with  the  abundance  of  “clean  water”  taxa  in  that 
year.  Also,  the  much  lower  levels  of  Chl -a  at  Cluny  in  1993  and  2006  as  compared  to  the 
1980s  supports  the  change  in  fauna  at  this  site  between  these  periods. 

A small  amount  of  epilithic  community  data  is  also  available.  In  this  case,  data  from 
October  1980  and  1981  (averaged)  were  compared  to  data  from  October  1994  and  2006. 
For  the  historical  data,  samples  from  Bowness  and  Bow  City  were  used  in  place  of 
Cochrane  and  Cluny.  Overall,  densities  were  lower  in  2006  than  for  the  other  two 
periods.  Carseland  had  the  highest  total  cell  densities  of  ail  sites  in  1994  and  2006,  while 
Bow  City  followed  by  Stier’s  Ranch  had  the  highest  densities  in  the  early  80s  (Figure 
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15).  Diatoms  made  up  the  greatest  percentage  of  density  at  most  sites  with  the  exception 
of  Bowness  in  80/81  and  Ronalane  in  1994  (Figure  16).  Green  algae  were  more 
prévalent  in  the  80/81  samples  than  in  1994  or  2006.  Overall,  cyanobacteria  contributed 
more  to  the  algal  density  at  Stier’s  Ranch  in  2006  than  historically. 
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4.0  GENERAL  DISCUSSION  AND  CONCLUSIONS 


Although  there  are  some  municipal  discharges  upstream  of  Cochrane  (i.e.,  Banff, 
Canmore)  the  sampling  site  at  Cochrane,  in  the  foothills  upstream  from  Calgary, 
expériences  low  nutrient  levels  and  showed  few  signs  of  nutrient  enrichment.  Epilithic 
algal  chlorophyll-rz  and  biomass  levels  were  low  and  a large  proportion  of  the  diatoms 
présent  were  taxa  which  prefer  low  levels  of  nutrients.  The  benthic  invertebrate 
community  at  this  site  was  impoverished  in  2006,  having  the  lowest  abundance  and 
diversity  of  ail  sites.  The  abundance  of  invertebrates  was  low,  in  particular  the  number 
and  variety  of  mayflies  and  stoneflies  was  much  lower  than  would  be  expected  at  a 
“clean  water”  site.  Rather  than  being  related  to  water  quality,  this  situation  is  likely  the 
resuit  of  rapidly  changing  flow  levels  from  the  operation  of  the  Ghost  Dam  for 
hydropower  génération.  Diurnal  fluctuations  in  discharges  from  the  dam  can  expose  and 
flood  the  riverbed  within  short  periods  of  time,  making  the  access  to  a constantly 
submerged  sampling  site  difficult.  In  addition,  and  perhaps  also  as  a resuit  of  flow 
régulation  ( e.g .,  Kirkwood  et  al.  2007),  Didymosphenia  geminata,  a colonial  diatom, 
frequently  forms  extensive  mats  of  cotton-like  material  at  this  site  and  contributes  further 
to  poor  habitat  for  invertebrates.  The  effects  of  these  perturbations  were  évident  in  some 
of  the  zoobenthic  samples  from  previous  years  (Anderson  1991).  Overall  benthic 
invertebrate  numbers  were  quite  variable  and  no  clear  temporal  trend  emerged. 

In  contrast,  Stier’s  Ranch,  immediately  downstream  from  two  wastewater  treatment 
plants,  showed  vivid  signs  of  organic  nutrient  enrichment.  It  had  the  highest  abundance 
of  organisms  and  the  richest,  most  diverse  benthic  community.  However,  the  numbers  of 
pollution-intolérant  taxa  tended  to  be  slightly  lower  at  this  site  than  at  sites  further 
downstream  from  Calgary’ s influence.  Nutrient  enrichment  was  also  apparent  in  the 
higher  epilithic  chlorophyll-a  levels,  larger  proportions  of  cyanobacteria,  and  higher 
number  of  diatom  taxa  preferring  greater  nutrient  concentrations  at  this  site.  The 
invertebrate  community  at  Stier’s  Ranch  in  2006  was  more  similar  to  samples  taken  in 
the  1980s  (i.e.,  after  implémentation  of  phosphorus  removal  at  Calgary’ s wastewater 
treatment  plants)  than  to  the  sample  taken  in  1993  (i.e.,  after  further  implémentation  of 
phosphorus  and  nitrogen  removal).  In  1993  taxonomie  diversity,  number  of  EPT  taxa, 
and  ratio  of  EPT  to  chironomids  were  higher,  and  the  biotic  index,  proportion  of 
nematodes  and  chironomids,  and  proportion  of  hydropsychids  to  Trichoptera  were  lower. 
These  are  ail  signs  that  conditions  were  better  in  1 993  than  in  the  other  years  that  were 
examined  in  detail.  While  invertebrates  clearly  respond  to  the  great  abundance  of  food, 
responses  to  diurnal  oxygen  déficit,  a recurring  summer  issue  in  the  Bow  River  (A. 
Sosiak,  2007  pers.  comm.),  were  not  apparent  at  the  erosional  sampling  location. 
However,  an  indication  of  such  response  to  stress  was  apparent  in  the  low  benthic 
invertebrate  diversity  indices  calculated  for  the  Stier’s  Ranch  depositional  site. 

Carseland  also  showed  the  effects  of  nutrient  enrichment.  Invertebrate  abundance, 
periphytic  algal  biomass,  cell  density,  and  chlorophyll-û  levels  were  relatively  high  at  this 
site  in  2006.  However,  the  numbers  of  pollution-sensitive  zoobenthic  taxa  were  slightly 
higher,  and  the  proportions  of  cyanobacteria  and  diatoms  preferring  high  nutrient  levels 
were  slightly  lower,  than  at  Stier’s  Ranch.  This  site  benefited  significantly  from  nutrient 
removal  in  the  1980s  (Sosiak  2002).  Based  on  available  data,  the  invertebrate 
community,  although  it  still  exhibits  responses  to  nutrient  enrichment,  has  remained  in 
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relatively  healthy  condition  since  that  time.  Given  its  location,  the  response  at  this  site  to 
Calgary’s  expansion  should  be  monitored  closely  in  the  future. 

While  invertebrate  abundance  was  lower  at  Cluny  than  at  Carseland  in  2006,  community 
composition  was  similar.  Epilithic  chlorophyll-Æ  was  also  lower  at  this  site  than  at 
Carseland;  in  the  1980s,  chlorophyll-a  levels  tended  to  be  higher.  Zoobenthic  abundance 
has  been  historically  higher  at  this  site;  however,  most  metrics  from  1985  suggest  that 
conditions  at  this  site  were  poor,  or  more  typical  of  depositional  conditions.  This  site 
exhibits  changes  over  time  which  suggest  a shift  from  a depositional  to  an  erosional 
environment.  Essentially,  zoobenthos  from  Cluny  has  become  more  like  that  from 
Carseland  and  less  like  that  from  Ronalane.  Such  changes  may  be  quite  independent  of 
water  quality  and  could  be  related  to  the  effect  of  flushing  flows  experienced  in  2005 
(Appendix  9)  and  to  différences  in  the  exact  sampling  location  historically  and  in  2006. 

In  2006,  Ronalane  was  defined  by  a large  proportion  of  burrowers  and  other  depositional 
fauna.  This  finding  is  consistent  with  this  site’ s location  near  the  mouth  of  the  Bow 
River  and  is  reflected  in  historical  conditions.  There  is  some  indication  that  this  site  may 
be  experiencing  détérioration.  While  flow  régulation  and  discharges  of  municipal 
effluents  hâve  well-documented  impacts  on  benthic  invertebrate  distribution,  the  effects 
of  other  potential  stressors,  such  as  irrigation  return  flows,  need  to  be  evaluated, 
particularly  at  lower  sites  on  the  Bow  River. 

Despite  the  temporal  changes  observed  at  some  of  the  sites,  the  consistent  way  in  which 
sites  are  separated  by  multivariate  analysis  suggests  that  many  of  the  basic  longitudinal 
différences  in  the  zoobenthic  community  of  the  Bow  River  remained  constant  over  the  23 
year  period  of  study.  This  is  an  indication  of  the  strong  influence  of  physical  and 
Chemical  ecoregional  différences  on  river  dynamics  and  benthic  invertebrate  distribution 
- an  observation  which  has  been  made  in  previous  studies  on  rivers  from  the  South 
Saskatchewan  River  basin  ( e.g .,  Culp  and  Davies  1982,  Anderson  1991,  Powell  2008). 

Aside  from  substrate  type,  flow  conditions  hâve  a definite  influence  on  benthic 
invertebrate  and  algal  population  density  and  taxa  distribution.  Biological  responses  to 
man-made  improvements  such  as  nutrient  load  réductions  from  municipal  wastewater 
discharges  may  be  overshadowed  by  responses  to  flow  fluctuations.  This  is  particularly 
so  for  the  Bow  River  benthic  invertebrate  data  set  which  had  too  few  collections  (time 
and  space)  to  adequately  tease  out  the  sometimes  confounding  influence  of  flow  and 
responses  to  municipal  wastewater  treatment  plant  improvements. 

In  order  to  build  on  these  data,  a regular  sampling  schedule  should  be  established.  More 
frequent  sampling  would  make  détermination  of  trends  in  biological  communities  due  to 
natural  and  human-induced  variations  in  flow  easier  to  détermine.  In  addition,  more 
intensive  spatial  sampling  would  facilitate  the  interprétation  of  longitudinal  trends  in 
water  quality  and  biota  as  they  relate  to  point-source  management  changes  in  the  basin. 

Overall,  benthic  invertebrate  and  epilithic  algae  data  collected  from  the  Bow  River 
provided  useful  and  complementary  information  about  longitudinal  and  temporal  trends 
at  Bow  River  long-term  monitoring  sites.  A large  number  of  metrics  was  used  and  it  is 
évident  that  they  each  provide  different  information,  some  being  more  indicative  of  the 
general  nutrient  status  of  a site,  others  of  the  physical  habitat  conditions.  These  different 
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but  important  features  need  to  be  recognized  and  metrics  for  analysis  selected 
accordingly. 

Detailed  identification  and  énumération  of  invertebrate  and  algal  samples  is  labour 
intensive,  but  essential  for  extracting  the  maximum  value  from  the  samples  collected. 
Rapid  assessment  techniques  generally  rely  on  invertebrate  data  at  higher  taxonomie 
levels  (i.e.,  family).  While  some  useful  information  can  be  derived  from  this  type  of 
analysis,  more  detailed  taxonomy  is  often  necessary  for  describing  certain  trends.  The 
need  for  greater  detail  was  apparent  for  habitat  preference  information.  The  évaluation  of 
metrics  based  on  other  classifications  such  as  behaviour  or  feeding  group  was  more 
successful  at  the  level  of  taxonomie  detail  available  for  this  study,  but  might  also  hâve 
been  improved  by  greater  taxonomie  detail.  Species-level  identification  proved 
important  for  evaluating  the  epilithic  algal  community  data,  as  this  allowed  the 
classification  of  diatoms  according  to  their  preferences  for  phosphorus  and  nitrogen 
levels. 
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An  Àssessmenî  ot  Benthic  Invertebrates  and  Epilithic  Algae  at  Long-term  Monitoring 

Sites  in  the  Bow  River  (Fali  2006) 


1.  Abondance 


2.  Richness 


Cochrane  Stier's  Carseland  Cluny  Ronalane 


Figure  2 Comparison  of  benthic  invertebrate  data  obtained  with  a Neill 
cylinder  and  kick  nets  (210  and  400pm  mesh  size)  in  erosional 
areas  from  the  Bow  River  (October  2006) 
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cylinder  and  kick  nets  (210  and  400pm  mesh  size)  in  erosional 
areas  from  the  Bow  River  (October  2006)  (con’t) 
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5.  Clinger  Species 


Figure  2 Comparison  of  benthic  invertebrate  data  obtained  with  a Neill 
cylinder  and  kick  nets  (210  and  400pm  mesh  size)  in  erosional 
areas  from  the  Bow  River  (October  2006)  (con’t) 
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Figure  2 Comparison  of  benthic  invertebrate  data  obtained  with  a Neil! 
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cylinder  and  kick  nets  (210  and  400jjm  mesh  size)  in  erosional 
areas  from  the  Bow  River  (October  2006)  (con’t) 
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Figure  2 Comparison  of  benthic  invertebrate  data  obtained  with  a Neill 
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areas  from  the  Bow  River  (October  2006)  (con’t) 
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13,  Contribution  of  Dominant  Order  (Insect) 
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cylinder  and  kick  nets  (210  and  400|jm  mesh  size)  in  erosional 
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Figure  2 Comparison  of  benthic  invertebrate  data  obtained  with  a Neill 
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19.  Biotic  Index  - Species 


Cochrane  Stier's  Carseland  Cluny  Ronalane 


Figure  2 Comparison  of  benthic  invertebrate  data  obtained  with  a Neill 
cylinder  and  kick  nets  (210  and  400pm  mesh  size)  in  erosional 
areas  from  the  Bow  River  (October  2006)  (con’t) 
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Figure  3 Bow  River  benthic  invertebrate  metrics  depicting  longitudinal 
trends  in  erosional  habitats  (October  2006) 
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5.  Shannon  Diversity  Index  2006 
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11.  Contribution  by  Dominant  Order  (Insect)  2006 
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Figure  3 Bow  River  benthic  invertebrate  metrics  depicting  longitudinal 
trends  in  erosional  habitats  (October  2006)  (con’t) 
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17.  Proportion  of  Predators  2006 
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22.  Burrower  Taxa  2006 
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Figure  3 Bow  River  benthic  invertebrate  metrics  depicting  longitudinal 
trends  in  erosional  habitats  (October  2006)  (con’t) 
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1.  ABUNDANCE 


Figure  4 Comparison  of  various  benthic  invertebrate  metrics  depicting 
longitudinal  trends  in  erosional  (E)  and  depositional  (D)  areas 
(October  2006) 
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4.  EPT  NUMBERS  AS  PERCENT  OF  TOTAL  NUMBER  OF  ORGANISMS 
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Figure  4 Comparison  of  various  benthic  Invertebrate  metrics  depictîng 
longitudinal  trends  in  erosional  (E)  and  depositional  (D)  areas 
(October  2006)  (con’t) 
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Figure  5 Comparison  of  habitat  preferences  of  invertebrates  in 

erosional  and  depositional  areas  from  the  Bow  River  (October 
2006) 
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Figure  6 Comparison  of  habits  of  invertebrates  in  erosional  and 
depositional  areas  from  the  Bow  River  (October  2006) 
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Figure  7 Plots  for  Discriminant  Functions  1, 2,  and  3 (each  Symbol  represents  data  from  1 Neill  sample 
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Figure  8 Summary  of  linear  régression  and  ANOVAs  on  discriminant  scores  by  year 
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Abundance  1993 


Figure  9-1  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006 

Total  refers  to  the  count  of  total  unique  taxa  per  site  (replicates  pooled) 
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Richness(Family)  1985 


Figure  9-2  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-3  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-4  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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EPT  Généra  2006 
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Figure  9-5  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-6  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-7  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-8  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985, 1993  and  2006  (con’t) 
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Figure  9-9  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-10  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-11  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Proportion  of  Oligochaetes  1985 


Proportion  of  Oligochaetes  1993 


Straight  line  indicates  the  mean  of  the  5 sites 

Figure  9-12  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-13  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-14  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-15  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-16  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985, 1993  and  2006  (con’t) 
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Figure  9-17  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-18  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-19  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-20  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-21  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-22  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  9-23  Comparison  of  benthic  invertebrate  metrics  at  erosional  sites  on  the 
Bow  River  for  1985,  1993  and  2006  (con’t) 
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Figure  10  Summary  of  metrics  for  epilithic  algae  collected  from  the  Bow  River 
in  2006 
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Figure  10  Summary  of  metrics  for  epilithic  algae  collected  from  the  Bow  River 
in  2006  (con’t) 
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Figure  11  Longitudinal  trends  in  epilithic  diatoms  classified  according  to  broad 
preference  ranges  for  phosphorus  and  nitrogen  concentrations 
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Figure  11  Longitudinal  trends  in  epilithic  diatoms  classified  according  to  broad 
preference  ranges  for  phosphorus  and  nitrogen  concentrations 

(con’t) 
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1.  Epilithic  Chlorophyll-a  (different  symbols  = different  sampler) 
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2.  Epilithic  Chlorophyll-a 


3.  Epilithic  Phaeophytin-a 


Figure  12  Longitudinal  trends  in  chlorophyll-a  in  epilithic  algal  samples  from 
the  Bow  River  (Oct  2006) 
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4.  Ratio  of  Epilithic  Chlorophyll-a  to  Phaeophytin-a 


Figure  12  Longitudinal  Trends  in  chlorophyll-a  in  epilithic  algal  samples  from 
the  Bow  River  (Oct  2006)  (con’t) 


An  Assessment  of  Benthic  Invertebrates  and  Epilithic  Algae  at  Long-term  Monitoring 
Sites  in  the  Bow  River  (Fall  2006) 


2.  Invertebrate  Shannon  Index  and  Epilithic  Chlorphyll-a  to 
Phaeophytin-a  Ratio 
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3.  Invertebrate  Richnessand  Epilithic  Algal  Richness 


Figure  13  Comparison  of  longitudinal  trends  for  a sélection  of  benthic 

invertebrate  and  epilithic  algal  metrics  in  the  Bow  River  (Oct  2006) 
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Figure  14  Comparison  of  longitudinal  trends  in  epilithic  chlorophyll-a  for  the 
Bow  River  in  1983-1986,  1993  and  2006 
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Figure  15  Comparison  of  epilithic  algal  density  in  the  Bow  River  in  1980-81, 
1994,  and  2006 
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Figure  16  Comparison  of  the  relative  density  (%)  of  epilithic  algal  taxa  in  the 
Bow  River  in  1980-81,  1994  and  2006 
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Appendix  2 Raw  benthic  invertebrate  data  from  the  Bow  River,  October  2006 
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Appendix  2 Raw  benthic  invertebrate  data  from  the  Bow  River,  October  2006  (con’t) 
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Appendix  3 Raw  epilithic  algal  data  for  the  Bow  River  (Fall  2006) 

Epilithic  algal  density  (number  of  units/cm2)  and  biomass  (micrograms/cm2)  at  sampling  sites  on  Bow  River 
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| Scenedesmus  opoliensis  P.  Richter  i 

| Scenedesmus  quadricauda  (Turpin)  Brebisson  ! 

[Scenedesmus  sempervirens  Chodat  j 

| Scenedesmus  sp  j 

a 
c n 

% 

I 

8 

.1 

55 

lCRYPTOPHYCEAE  j 

| Cryptomonas  sp  ; 

| Rhodomonas  minuta  Skuja  i 

IEUGLENOPHYCEAE  | 

|Eugf/ena  cf.  minuta  Prescott  j 
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Appendix  4 Metrics  derived  from  benthic  invertebrate  data  collected  with  a Neill 
cylinder  and  kick  nets  (210  and  400  pm  mesh  size)  in  erosional  and 
depositional  areas  from  the  Bow  River  in  Fall  2006 


Cochrane 

Stier's 

Carseland 

Metrics/lndices 

NEILL 

210E 

400E 

210D 

NEILL 

210E 

400E 

21 0D 

NEILL 

210E 

400E 

210D 

Abundance  (total  count) 

667 

2112 

337 

1973 

14811 

36240 

18438 

13940 

8623 

73675 

16571 

61634 

Richness  (total  species) 

43 

39 

39 

43 

70 

45 

45 

48 

75 

68 

56 

59 

Total  Taxa  (qenus  and  hiqher) 

41 

35 

36 

39 

68 

45 

45 

47 

74 

66 

56 

51 

Total  Taxa  (family  and  higher) 

25 

24 

26 

27 

37 

29 

25 

28 

39 

36 

28 

29  ! 

Simpson's  Diversity  Index  D 

0.80 

0.90 

0.90 

0.88 

0.92 

0.88 

0.89 

0.78 

0.90 

0.77 

0.86 

0.37 

Shannon  Diversity  Index  H' 

2.18 

2.60 

2.78 

2.57 

2.84 

2.55 

2.64 

2.16 

2.78 

2.20 

2.59 

1.05 

Number  of  Coleoptera  généra 

2 

2 

2 

2 

4 

2 

3 

1 

3 

3 

4 

3 

Number  of  Diptera  généra  (non 
Chironomidae) 

3 

2 

4 

! 

2 

2 

2 

1 

4 

4 

3 

1 

Number  of  Chironomid  généra 

14 

10 

10 

11 

22 

15 

16 

17 

27 

24 

21 

17 

Number  of  Hemiptera  qenera 

2 

3 

2 

3 

3 

2 

1 

2 

1 

2 

0 

5 

Number  of  Odonata  généra 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

Number  of  Ephemeroptera  qenera 

3 

1 

2 

2 

6 

4 

7 

5 

8 

7 

9 

4 

Number  of  piecoptera  généra 

0 

0 

1 

0 

2 

1 

1 

1 

4 

4 

3 

0 

Number  of  T richoptera  qenera 

3 

1 

31 

2 

8 

6 

5 

3 

9 

6 

7 

1 

Number  of  EPT  généra 

6 

2 

6 

4 

16 

11 

13 

9 

21 

17 

19 

5 

Number  of  Coleoptera  families 

2 

2 

2 

2 

2 

2 

2 

1 

2 

3 

2 

2 

Number  of  Diptera  (non 
Chironomidae)  families 

2 

2 

4 

1 

2 

2 

2 

1 

3 

4 

3 

1 

Number  of  Chironomid  subfamilies 

5 

5 

5 

6 

5 

5 

5 

5 

5 

5 

5 

5 

Number  of  Hemiptera  families 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

1 

Number  of  Odonata  families 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

Number  of  Ephemeroptera  families 

2 

1 

2 

2 

4 

4 

4 

4 

6 

5 

6 

3 

Number  of  Piecoptera  families 

0 

0 

1 

0 

2 

1 

1 

1 

2 

2 

2 

0 

Number  of  T richoptera  families 

3 

1 

3 

2 

6 

5 

4 

2 

7 

4 

5 

1 

Number  of  EPT  families 

5 

2 

6 

4 

12 

10 

9 

7 

15 

11 

13 

4 

Number  of  ETO  families 

5 

2 

5 

4 

10 

9 

8 

6 

13 

9 

12 

4 

Proportion  of  Oliqochaetes 

0.2908546 

0.2386364 

0.379822 

0.4460213 

0.093536 

0.013797 

0.027064 

0.457102 

0.1844455 

0.4140618 

0.1705992 

0.7957459 

Proportion  of  Gastropods 

0.002099 

0 

0.0029674 

0.0015205 

0.033839 

0.017384 

0.056676 

0.04462 

0.0020643 

0.0021853 

0.0024742 

0.0001785 

Proportion  of  Pelecypods 

0 

0.0004735 

0 

0 

0.00027 

2.76E-05 

0.003525 

0.000574 

4.639E-05 

0 

0 

0 

Proportion  of  Amphipods 

0 

0 

0 

0.0010137 

0.0005 

5.52E-05 

0.001193 

0.000646 

2.319E-05 

0.0005022 

0 

0.0009086 

Proportion  of  Coleoptera 

0.0083958 

0.0137311 

0.0356083 

0.0136847 

0.050988 

0.057892 

0.058412 

7.17E-05 

0.0276019 

0.0113607 

0.0260696 

0.0009086 

Proportion  of  Diptera 

0.0722639 

0.0677083 

0.1869436 

0.1282311 

0.466195 

0.570088 

0.409914 

0.267791 

0.6079837 

0.474245 

0.597912 

0.0916215 

Proportion  of  Diptera  (non 
Chironomidae) 

0.002099 

0.0047348 

0.0207715 

0.0005068 

0.120326 

0.064018 

0.214665 

0.002941 

0.000719 

0.003149 

0.0126124 

6.49E-05 

Proportion  of  Chironomids 

0.0701649 

0.0629735 

0.1661721 

0.1277243 

0.345869 

0.506071 

0.195249 

0.264849 

0.6072646 

0.471096 

0.5852996 

0.0915566 

Proportion  of  Hemiptera 

0.0245877 

0.1917614 

0.1958457 

0.1084643 

0.009344 

0.00091 1 

0.077069 

0.000359 

4.639E-05 

0.0001222 

0 

0.0058085 

Proportion  of  Odonata 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6.035E-05 

0 

Proportion  of  Ephemeroptera 

0.0041979 

0.0023674 

0.0652819 

0.0010137 

0.038659 

0.07489 

0.166395 

0.013199 

0.0310115 

0.0390635 

0.0698208 

0.0008112 

Proportion  of  Piecoptera 

0 

0 

0.0029674 

0 

6.75E-05 

0.000276 

0.000542 

7.17E-05 

0.0005567 

0.0008687 

0.007483 

0 

Proportion  of  Trichoptera 

0.0026987 

0.0056818 

0.0237389 

0.003041 1 

0.131129 

0.239294 

0.186137 

0.002869 

0.065943 

0.0185816 

0.1085632 

8.1 12E-05 

Proportion  of  EPT 

0.0068966 

0.0080492 

0.0919881 

0.0040547 

0.169856 

0.314459 

0.353075 

0.016141 

0.0975112 

0.0585137 

0.1858669 

0.0008924 

Ratio  of  EPT  to  Chironomidae 

0.0982906 

0.1278195 

0.5535714 

0.031746 

0.491099 

0.621374 

1 .808333 

0.060943 

0.1605745 

0.1242077 

0.3175585 

0.0097466 

contribution  by  dominant 
family(insects  only) 

0.644385 

0.6818182 

0.3837209 

0.5039841 

0.669452 

0.604323 

0.456236 

0.941726 

0.8292837 

0.8713869 

0.738246 

0.9233159 

Proportion  of  Filterers 

0.0335832 

0.0648674 

0.0534125 

0.0968069 

0.298594 

0.309161 

0.444896 

0.128336 

0.2307889 

0.1281303 

0.206264 

0.0442937 

FiltererTaxa 

7 

9 

5 

8 

13 

11 

10 

11 

12 

10 

10 

9 

Proportion  of  Gatherers 

0.4818591 

0.5127841 

0.5727003 

0.6482514 

0.297339 

0.323538 

0.244224 

0.650287 

0.5330411 

0.7323108 

0.640939 

0.9004932 

Gatherer  Taxa 

21 

14 

11 

14 

23 

15 

14 

20 

29 

27 

20 

23 

Proportion  of  Parasites 

0.391904 

0.1690341 

0.0445104 

0.0983274 

0.14262 

0.001104 

0.001085 

0.00581 1 

0.0451372 

0.0077367 

0.0067588 

0.0188208 

Parasite  Taxa 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Proportion  of  Piercing  Herbivores 

0.0002999 

0.0004735 

0.0029674 

0.0187532 

4.05E-05 

0 

0 

7.17E-05 

0 

6.787E-05 

0 

0.0005192 

Piercing  Herbivore  Taxa 

1 

1 

1 

3 

1 

0 

0 

1 

0 

1 

0 

2 

Proportion  of  Predators 

0.0734633 

0.2466856 

0.2908012 

0.1322859 

0.107849 

0.187086 

0.148877 

0.052869 

0.0529771 

0.0326026 

0.0519582 

0.0045267 

PredatorTaxa 

9 

12 

12 

11 

21 

11 

10 

10 

19 

15 

13 

17 

Proportion  of  Scrapers 

0.0047976 

0.0061553 

0.0237389 

0.0050684 

0.108038 

0.086065 

0.14405 

0.159613 

0.0361376 

0.0330641 

0.0525617 

0.0247266 

ScraperTaxa 

3 

2 

4 

5 

7 

4 

8 

4 

8 

9 

7 

6 

Proportion  of  Shredders 

0.014093 

0 

0.0118694 

0.0005068 

0.045519 

0.093046 

0.016867 

0.003013 

0.1019182 

0.0660875 

0.0415183 

0.0066197 

Shredder  T axa 

3 

0 

4 

1 

2 

2 

2 

2 

4 

4 

4 

2 

Scraper/Filterer  Ratio 

0.1428571 

0.0948905 

0.4444444 

0.052356 

0.361823 

0.278383 

0.323784 

1.243712 

0.1565829 

0.2580508 

0.2548274 

0.5582418 

Shredders/Total  Ratio 

0.014093 

0 

0.0118694 

0.0005068 

0.045519 

0.093046 

0.016867 

0.003013 

0.1019182 

0.0660875 

0.0415183 

0.0066197 

Total  Insects 

75 

594 

172 

502 

10314 

34187 

16566 

3964 

6322 

40097 

13421 

6116 

Total  Insect  Taxa 

29 

23 

27 

25 

49 

32 

35 

31 

57 

52 

48 

39 

Proportion  of  Burrowers 

0.0187166 

0 

0.005814 

0.1195219 

0.023404 

0.00819 

0.000302 

0.204844 

0.0926981 

0.041649 

0.0730944 

0.0168411 

BurrowerTaxa 

3 

0 

1 

1 

6 

2 

1 

3 

8 

6 

6 

3 

Proportion  of  Climbers 

0.1042781 

0.0420875 

0.0348837 

0.249004 

0.01687 

0.003861 

0.008874 

0.067104 

0.112345 

0.1498865 

0.0502943 

0.3019948 

Climber  Taxa 

4 

3 

2 

4 

4 

3 

4 

3 

4 

4 

4 

3 

Proportion  of  Clingers 

0.1657754 

0.0488215 

0.1802326 

0.069721 1 

0.618553 

0.604411 

0.671013 

0.509586 

0.3892053 

0.2575754 

0.3457268 

0.3616743 

ClingerTaxa 

7 

5 

8 

6 

21 

15 

16 

11 

24 

22 

21 

10 

Proportion  of  Sprawlers 

0.3903743 

0.1801347 

0.2848837 

0.0836653 

0.296149 

0.340363 

0.111433 

0.216196 

0.367407 

0.4859466 

0.4717234 

0.2532701 

Sprawler  Taxa 

9 

8 

9 

6 

8 

7 

8 

8 

15 

13 

11 

6 

Proportion  of  Swimmers 

0.3208556 

0.7289562 

0.494186 

0.4780876 

0.045024 

0.043174 

0.208379 

0.00227 

0.0383447 

0.0649425 

0.059161 

0.0662198 

SwimmerTaxa 

6 

7 

7 

8 

10 

5 

6 

6 

6 

7 

6 

17 

Proportion  with  Depositional  pref. 

0.0935829 

0.3265993 

0.0988372 

0.6733068 

0.088885 

0.010472 

0.046239 

0.256307 

0.1142749 

0.0842208 

0.0530512 

0.2671681 
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Appendix  4 Metrics  derived  from  benthic  invertebrate  data  collected  with  a Neill 
cylinder  and  kick  nets  (210  and  400  pm  mesh  size)  in  erosionai  and 
depositionai  areas  from  the  Bow  River  in  Fall  2006  (con’t) 


Cochrane 

Stier's 

Carseland 

Metrics/lndices 

NEILL 

210E 

400E 

210D 

NEILL 

210E 

400E 

210D 

NEILL 

21 0E 

400E 

210D 

Depositionai  Taxa 

5 

9 

7 

10 

11 

6 

6 

9 

12 

13 

9 

15 

Proportion  with  Erosionai  prêt. 

0.4304813 

0.1767677 

0.2906977 

0.1314741 

0.657547 

0.715155 

0.579983 

0.270434 

0.5507466 

0.6042846 

0.6688026 

0.3629823 

Erosionai  Taxa 

9 

8 

8 

7 

20 

14 

15 

10 

25 

22 

20 

8 

Proportion  Erosionai  - Depositionai 
prêt. 

0.1631016 

0.0774411 

0.2325581 

0.0717131 

0.124738 

0.1304 

0.235784 

0.046418 

0.0844406 

0.0950196 

0.117726 

0.0179856 

Erosionai  - Depositionai  Taxa 

11 

4 

8 

7 

9 

7 

9 

6 

12 

10 

13 

6 

Proportion  with  Lentic  prêt. 

0.3101604 

0.4074074 

0.3488372 

0.123506 

0.108101 

0.104192 

0.128214 

0.416751 

0.2023538 

0.1797641 

0.1007376 

0.3453237 

LenticTaxa 

3 

1 

3 

1 

5 

3 

4 

5 

3 

4 

3 

9 

Proportion  with  Lotie  (general)  prêt. 

0.0026738 

0.0117845 

0.0290698 

0 

0.020728 

0.039781 

0.009779 

0.010091 

0.048184 

0.036711 

0.0596826 

0.0065402 

Lotie  (general)  Taxa 

1 

1 

1 

0 

4 

2 

1 

1 

5 

3 

3 

1 

Proportion  with  Margin  habitat  prêt. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Margin  Taxa 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Biotic  Index  (species) 

7.81 

7.09 

6.21 

7.86 

5.84 

5.16 

5.10 

7.41 

6.23 

6.78 

5.91 

7.81 

Bl  Arthropods  Only 

6.88 

6.48 

4.87 

6.47 

5.23 

5.05 

4.77 

7.19 

5.80 

5.93 

5.47 

7.06 

Bl  Insects  Only 

5.28 

5.01 

4.75 

5.50 

5.12 

4.99 

4.73 

6.84 

5.73 

5.86 

5.47 

6.31 

Family  Biotic  Index  (Modified) 

7.16 

6.35 

6.15 

7.11 

5.95 

5.38 

5.19 

7.28 

6.54 

6.89 

6.03 

7.81 

FBI  Arthropods  only 

5.59 

5.09 

4.93 

5.76 

5.62 

5.32 

4.99 

7.07 

6.20 

6.10 

5.62 

6.57 

FBI  Insects  Only 

5.78 

5.18 

4.95 

5.80 

5.63 

5.33 

4.99 

7.10 

6.21 

6.11 

5.62 

6.56 

Number  of  Intolérant  Taxa  (0,1,2) 

5 

3 

4 

2 

9 

5 

7 

4 

12 

10 

11 

2 

BMW  P Score  (Sensitivity) 

90 

68 

107 

85 

162 

119 

115 

110 

180 

160 

150 

78 

Average  Score  Per  T axon  (ASPT) 

5.00 

4.25 

5.10 

4.47 

4.91 

5.17 

5.00 

4.78 

5.29 

5.33 

5.77 

4.11 
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Appendix  4 Metrics  derived  from  benthic  invertebrate  data  collected  with  a Neill 
cylinder  and  kick  nets  (210  and  400  pm  mesh  size)  in  erosional  and 
depositional  areas  from  the  Bow  River  in  Fall  2006  (con’t) 


Cluny 

Ronalane 

Metrics/lndices 

NEILL 

21 0E 

400E 

210D 

NEILL 

21 0E 

400E 

21 0D 

Abundance  (total  count) 

4396 

23542 

9027 

37274 

7447 

34788 

26915 

24881 

Richness  (total  species) 

65 

56 

60 

39 

68 

71 

71 

46 

Total  Taxa  (qenus  and  higher) 

65 

56 

60 

34 

67 

70 

71 

39 

Total  Taxa  (family  and  higher) 

37 

32 

36 

18 

43 

40 

41 

27 

Simpson's  Diversity  Index  D 

0.90 

0.70 

0.87 

0.27 

0.84 

0.80 

0.77 

0.71 

Shannon  Diversity  Index  H' 

2.71 

1.93 

2.59 

0.63 

2.35 

2.13 

2.01 

1.69 

Number  of  Coleoptera  généra 

4 

2 

2 

1 

3 

4 

3 

2 

Number  of  Diptera  généra  (non 
Chironomidae) 

5 

2 

3 

1 

5 

8 

10 

! 3 

Number  of  Chironomid  généra 

19 

17 

16 

12 

18 

20 

17 

7 

Number  of  Hemiptera  généra 

2 

0 

2 

5 

2 

1 

2 

6 

Number  of  Odonata  généra 

0 

0 

0 

0 

1 

0 

1 

1 

Number  of  Ephemeroptera  généra 

9 

8 

9 

3 

5 

9 

7 

2 

Number  of  Plecoptera  généra 

3 

3 

4 

Ô 

0 

1 

2 

0 

Number  of  T richoptera  généra 

9 

8 

8 

2 

6 

5 

8 

1 

Number  of  EPT  généra 

21 

19 

21 

5 

11 

15 

17 

3 

Number  of  Coleoptera  families 

2 

2 

2 

1 

2 

3 

2 

2 

Number  of  Diptera  (non 
Chironomidae)  families 

4 

2 

3 

1 

4 

5 

6 

2 

Number  of  Chironomid  subfamilies 

5 

5 

4 

4 

4 

4 

5 

4 

Number  of  Hemiptera  families 

1 

0 

1 

1 

1 

1 

1 

1 

Number  of  Odonata  families 

0 

0 

0 

0 

1 

0 

1 

1 

Number  of  Ephemeroptera  families 

6 

4 

6 

2 

4 

5 

4 

2 

Number  of  Plecoptera  families 

2 

2 

2 

0 

0 

1 

1 

0 

Number  of  Trichoptera  families 

7 

6 

6 

2 

4 

3 

5 

1 

Number  of  EPT  families 

15 

12 

14 

4 

8 

9 

10 

3 

Number  of  ETO  families 

13 

10 

12 

4 

9 

8 

10 

4 

Proportion  of  Oligochaetes 

0.236113 

0.522817 

0.143833 

0.851666 

0.3183393 

0.3382712 

0.444654 

0.3125678 

Proportion  of  Gastropods 

0 

0 

3.69E-05 

0 

8.057E-05 

5.749E-05 

0.001214 

0 

Proportion  of  Pelecypods 

4.55E-05 

0.000581 

7.39E-05 

0 

0.0016113 

0.0032962 

0.006626 

0.0107713 

Proportion  of  Amphipods 

0 

0.000566 

3.69E-05 

0 

0.00333 
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Appendix  4 Metrics  derived  from  benthic  invertebrate  data  collected  with  a Neill 
cylinder  and  kick  nets  (210  and  400  pm  mesh  size)  in  erosional  and 
depositional  areas  from  the  Bow  River  in  Fall  2006  (con’t) 
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Appendix  5 Benthic  invertebrate  data  collected  with  a Neill  cylinder  sampler  from  the  Bow  River  (1983  to  2006) 
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Appendix  5 Benthic  invertebrate  data  collected  with  a Neill  cylinder  sampler  from  the  Bow  River  (1983  to  2006) 
(con’t) 
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Appendix  6 Summary  of  results  of  the  Bootstrap  assessment  of  différences  among  sampling  approaches 
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Appendix  7 Summary  of  ANOVAS  performed  on  2006  benthic  invertebrate 
data  to  assess  longitudinal  trends  in  the  Bow  River 


Test  p for  ANOVA  Bonferroni  (p<0.05) 

Metric  Tested 

Cochrane 
Stiers  Ranch 
Carseland 
Cluny 
Ronalane 

Abundance 

0.000 

(i) 

Richness 

0.000 

Simpson  Diversity  Index 

0.017 

Shannon  Diversity  Index 

0.000 

EPT  généra 

0.000 

EPT  to  chironomid  ratio 

0.000 

Contribution  by  Dominant  order  (insect) 

0.000 

Number  of  Filter  Taxa 

Ô.ÔOÔ 

Proportion  of  Filterers 

0.000 

Proportion  of  Gatherers 

! 0.000 

Proportion  of  Shredders 

0.427 

no  sig 

inifice 

rnt  diff 

erenc 

;es 

Proportion  Scrapers 

0.000 

Proportion  Predators 

0.000 

Proporation  Swimmers 

0.000 

Number  of  Clinger  Taxa 

0.000 
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Appendix  7 Summary  of  ANOVAS  performed  on  2006  Benthic  Invertebrate 
Data  to  Assess  Longitudinal  Trends  in  the  Bow  River(con’t) 


Test  P for  ANOVA  Bonferroni  (p<0.05) 

Metric  Tested 

Cochrane 
Stiers  Ranch 
Carseland 
Cluny 
Ronalane 

Number  of  Clinger  Taxa 

0.000 

Proportion  Sprawlers 

0.005 

Number  of  BurrowerTaxa 

0.000 

Proportion  Burrowers 

0.000 

Proportion  with  Depositional  Preference  (Historié  Data) 

0.000 

Number  of  Pollution  Intolérant  Taxa 

o.ôôô 

Family  Biotic  Index 

0.000 

Sensitivity  (BWMP)  Score 

O.OOO 

Average  Score  per  Taxon 

0.000 

( 1 ) Horizontal  Unes  join  2 sites  that  differ  significantly  (p<0. 05) 
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Appendix  8 Summary  of  results  of  Step-wise  Discriminant  Analysis  on 
historical  benthic  invertebrate  Neill  cylinder  data  for  the 
Bow  River 


Function 

Eigenvalue 

% of 
Variance 

Cumulative 

% 

Canonical 

Corrélation 

1 

12.934 

43.4 

43.4 

0.963 

2 

8.872 

29.8 

73.2 

0.948 

3 

6.845 

23.0 

96.2 

0.934 

! 4 

1.142 

3.8 

100.0 

0.730 

Corrélations  for  individual  taxa  bv  function 


Taxon 

Function  -1 

Function  -2 

Function  -3 

Function  -4 

Chironomini 

0.33 

0.001 

0.043 

-0.011 

Tubificidae 

0.248 

-0.04 

-0.024 

0.109 

Simulidae 

0.203 

-0.108 

-0.175 

-0.062 

Coleoptera 

0.185 

-0.024 

0.098 

-0.154 

Tanypodinae 

0.163 

-0.001 

0.009 

0.115 

Hydropsychidae 

0.121 

-0.09 

-0.005 

-0.005 

Stenonema 

0.122 

0.271 

0.191 

0.115 

Hydra 

-0.208 

-0.255 

-0.053 

0.022 

T ricorythodes 

0.211 

0.218 

0.136 

-0.027 

Ephemerellidae 

0.068 

-0.162 

-0.027 

-0.117 

Orthocladiinae 

0.078 

-0.091 

-0.013 

-0.063 

Helicopsychidae 

0.192 

-0.201 

0.36 

0.23 

Hydracarina 

-0.026 

0.067 

-0.088 

-0.006 

Tanytarsini 

0.038 

0.06 

0.077 

-0.005 

Naididae 

-0.074 

-0.135 

-0.013 

-0.337 

Hydroptilidae 

-0.042 

0.099 

0.163 

0.329 

Hydropsyche 

0.187 

0.235 

-0.239 

-0.324 

Chironomidae 

0.014 

0.003 

0.123 

0.259 

Psychomyia 

0.226 

-0.175 

0.127 

0.23 

NEMATODA 

0.045 

0.002 

0.124 

-0.22 

Baetidae 

0.112 

-0.061 

-0.03 

-0.208 

OSTRACODA 

0.024 

-0.007 

-0.124 

0.188 

COPEPODA 

0.022 

-0.057 

-0.15 

0.165 

Enchytraeidae 

-0.124 

-0.036 

-0.001 

-0.164 

Variables  ordered  by  absolute  size  of  corrélation  within  function.  

Highlighted  numbers  indicate  the  largest  absolute  corrélation  between  each  variable 
and  any  discriminant  function 
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Appendix  9 Discussion  on  the  influence  of  discharge  on  variability  in 
the  historical  benthic  invertebrate  data  from  the  Bow  River 

CONSIDERATIONS  ABOUT  THE  INFLUENCE  OF  RIVER  FLOW  ON  BENTHIC 
INVERTEBRATE  DISTRIBUTION 

Différences  in  benthic  fauna  among  years  may  be  the  resuit  of  changes  to  water  quality 
over  time.  However,  the  effect  of  flow  on  invertebrate  community  structure  must  be 
considered  as  well.  Hydrometric  data  which  correspond  to  the  period  studied  are 
available  from  the  Water  Survey  of  Canada  on-line  archives  for  a limited  number  of  Bow 
River  sites.  Bow  River  below  Bearspaw  Dam  (05BH009),  downstream  of  Carseland 
Dam  (05BM002),  and  near  the  mouth  (at  the  confluence  of  the  Oldman  River) 

(05BN012)  were  chosen  to  represent  flow  conditions  over  the  range  of  sites  where 
invertebrates  were  collected.  These  stations  roughly  correspond  to  the  Cochrane, 
Carseland,  and  Ronalane  sampling  locations. 

Flows  from  the  dates  sampled  ranged  from  45  to  103  m3/s  downstream  from  the 
Bearspaw  Dam.  Downstream  from  the  Carseland  Dam,  daily  flows  on  the  sample  dates 
ranged  more  widely,  between  41  and  155  m3/s.  Near  the  confluence  with  the  Oldman 
River,  samples  were  collected  over  the  greatest  range  of  flows,  between  13  and  143  m3/s. 
In  ail  cases,  the  highest  flows  were  in  September  of  1993,  while  the  date  of  the  lowest 
flow  was  different  for  ail  sites  (Figure  A9-1).  The  highest  daily  flow  for  each  site 
occurred  in  the  sample  month  with  the  highest  mean  monthly  flows.  However,  the  lowest 
mean  monthly  flows  did  not  correspond  to  lowest  daily  flows  (Figure  A9-2). 

Hydrometric  conditions  prior  to  the  day  or  month  in  which  the  invertebrates  were 
sampled  can  hâve  an  important  effect  on  the  benthic  invertebrate  community.  The 
scouring  action  of  high  flows  may  reduce  invertebrate  abundance  and  diversity  in  the 
short-term,  but  by  removing  algal  over-growth  and  sédiments  may  leave  a substrate  that 
is  more  favourable  for  colonization  by  clinging  and  sprawling  taxa  that  are  associated 
with  “cleaner”  water  conditions.  Daily  flows  for  April  through  October  for  the  years 
sampled  are  graphed  in  the  Appendix  (Figures  A 9-3).  These  data  were  examined  for 
incidents  of  “high”  flow,  arbitrarily  defined  as  flows  over  300  m3/s.  The  first  such  peak 
occurred  downstream  of  the  Carseland  Dam  and  near  the  confluence  with  the  Oldman 
River  in  mid-September  of  1985,  just  two  weeks  before  invertebrate  sampling  occurred. 
On  this  occasion,  flows  rose  from  40  m3/s  to  about  330  m3/s  over  one  day.  In  early  June 
of  1986,  spring  flows  peaked  at  over  400  m3/s  at  ail  three  hydrometric  stations.  Flows  in 
1993  were  generally  high;  measurements  over  300  m3/s  occurred  numerous  times  at  the 
two  downstream  stations  throughout  June,  July  and  August.  The  maximum  daily  flow 
downstream  of  the  Carseland  Dam  was  490  m3/s  during  this  period,  while  flows  near  the 
mouth  peaked  at  477  m3/s.  The  highest  levels  of  ail  the  years  sampled  occurred  in  2006, 
when,  in  mid-June,  maximum  flows  reached  820  m3/s  downstream  of  the  Carseland  Dam 
and  778  m3/s  near  the  confluence. 

Although  no  benthic  samples  were  collected  in  2005,  it  is  important  to  note  that  in  June 
of  that  year,  the  Bow  River  experienced  a significant  flood  event.  No  flows  are  available 
from  the  station  below  the  Bearspaw  Dam,  but  the  highest  daily  flow  downstream  at 
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Calgary  (05BH004)  was  602  m3/s.  Flows  were  highest  downstream  of  the  Carseland 
Dam,  reaching  1710  m3/s.  Near  the  mouth  of  the  Bow  River,  the  maximum  daily  flow 
was  1530  m3/s.  Flows  of  these  magnitudes  would  be  expected  to  cause  physical 
alteration  to  the  bed  and  banks  of  the  river,  décimation  of  macrophyte  and  algae 
populations,  and  significant  drift  and  displacement  of  benthic  invertebrates. 

Ecological  conditions  in  the  Bow  River  in  the  fall  of  2004  and  2005  were  compared  in  a 
project  by  University  of  Calgary  Environmental  Sciences  students  (ENSC502  2006). 
From  an  examination  of  four  sites  within  the  City  of  Calgary,  they  concluded  that  while 
epilithic  algae  had  already  begun  to  recover  to  pre-flood  levels,  aquatic  macrophyte 
biomass,  sédiment  nutrient  levels,  and  benthic  invertebrate  diversity  (richness)  were  ail 
lower  than  in  the  previous  year.  Given  the  relatively  short  life-spans  of  most  benthic 
invertebrates,  recovery  should  be  expected  by  the  following  year  (i.e.  by  the  2006 
sampling  dates).  Indeed,  anecdotal  reports  from  fly-fishermen  claim  that  the  trout  fishery 
between  Calgary  and  Carseland  was  excellent  in  2006  (Big  Rock  Fly  Guides  2007), 
suggesting  an  adéquate  supply  of  invertebrate  prey. 
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Daily  Flow  d/s  Bearspaw  Dam  for  Sample  Dates 


1983  1984  1985  1986  1987  1989  1993  2006 


Daily  Flow  d/sCarseland  Dam  for  Sample  Dates 


1983  1984  1985  1986  1987  1989  1993  2006 


Line  indicates  médian 


Figure  A 9-1  Daily  flows  in  the  Bow  River  which  correspond  to  benthic 
invertebrate  sampling  dates,  1983-2006 
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Mean  Monthly  Flow  d/s  Bearspaw  Dam  for  Sample  Months 


Mean  Monthly  Flow  d/sCarseland  Dam  for  Sample  Months 


Mean  Monthly  Flow  Near  the  Mouth  for  Sample  Months 


Line  indicates  médian 

Figure  A 9-2  Monthly  flows  in  the  Bow  River  which  correspond  to 
benthic  invertebrate  sampling  months,  1983-2006 
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S/CUJ  D S/Eiu  o s/cui  o 


Below  Bearspaw  1983 


Confluence  1983 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River 
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s/Eiu  O S/Eui  o S/£UI  o 


Below  Bearspaw  1984 


DS  Carseland  1984 


Confluence  1984 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 
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s/eui  0 S/Eui  o S/Cui  o 


Below  Bearspaw  1985 


DS  Carseland  1985 


Confluence  1985 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 
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Below  Bearspaw  1986 


DS  Caresland  1986 


Confluence  1986 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 


An  Assessment  of  Benthic  Invertebrates  and  Epilithic  Algae  at  Long-term  Monitoring  153 

Sites  in  the  Bow  River  (Fall  2006) 


s/eui  o s/eui  o s/eiu  o 


Below  Bearspaw  1987 


DS  Carseland  1987 


Confluence  1987 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 
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s/eui  0 s/eui  o S/Gi u o 


Below  Bearspaw  1989 


DS  Carseland  1989 


Confluence  1989 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 
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Below  Bearspaw  1993 


DS  Caresland  1993 


Confluence  1993 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 
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s/eui  O S/Eui  o s/£iu  o 


Below  Bearspaw  2006 


DS  Ca resla nd  2006 


Confluence  2006 


Pink  square  indicates  benthic  invertebrate  sampling  day 

Figure  A 9-3  Daily  Flows  at  several  locations  in  the  Bow  River  (con’t) 
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AppendixIO  Comparisons  among  years:  examples  of  results  of 

Bonferroni  test  performed  on  samples  collected  in  1985, 
1993,  and  2006 

Example  applies  to  selected  metrics  from  Stier's  Ranch 

Metric  Tested  1985  1993  2006 

Abundance  

Richness 


Simpson  Diversity  Index 
Shannon  Diversity  Index 
No.  Predatory  Taxa 
No.  Scraper  Taxa 


No.  Shredder  taxa 
No.  Filter  Taxa 


No.  Gatherer  Taxa 


No.  Clinger  Taxa 


No.  Climbing  Taxa 


No.  Burrowing  Taxa 


No.  Sprawling  Taxa 


No.  Swimming  Taxa 


No.  Erosional  Taxa 


No.  Lentic  Taxa 


No.  Intolérant  Taxa 


Total  No.  Families 


No.  Ephemeroptera  Généra 


No.  Plecoptera  Généra 


No.  Trichoptera  Généra 


Full  line  joins  years  significantly  (p<0.05)  different 
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